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Abstract 



The consequences of application of new set of criteria, proposed in our previous works, for 
the improvement of a jet energy calibration accuracy with the process ll pp — ► 7 + jet + X" at 
Tevatron and for a reduction of background events contribution are studied. The efficiencies of 
the used selection criteria are estimated. The distributions of these events over and r^ 6 * are 
presented. The features of "7 + jet" events in the central calorimeter region of the DO detector 
(| 77 1 <0.7) are investigated. 

It is also shown that the samples of "7 + jet" events, selected with the cuts used for the 
jet energy calibration, may have the statistics sufficient for determining the gluon distribution 
function of a proton in the region of 2 x 1CT 3 < x < 1.0 and the values of Q 2 by one order higher 
than that reached in the experiments at HERA. 

Monte Carlo events produced by the PYTHIA 5.7 generator are used for this aim. 



1. INTRODUCTION. 



Setting an absolute energy scale for a jet, detected mostly by hadronic and electromagnetic 
calorimeters (HCAL and ECAL), is an important task for any pp or pp collider experiment (see 
e.g. [1-8]). 

The main goal of this work is to find out the selection criteria for "pp — - > 7 + jet + X" 
events (we shall use in what follows the abbreviation "7 + jet" for them) that would lead to 
the most precise determination of the transverse momentum of a jet (i.e. P^ et ) via assigning a 
photon P{' to a signal produced by a jet. Our study is based on the "7 + jet" events generated 
by using PYTHIA 5.7 J9|. Their analysis was done on the "particle level" (in the terminology 
of [1]), i.e. without inclusion of detector effects. The information provided by this generator 
is analyzed to track starting from the parton level (where parton-photon balance is supposed to 
take place in a case of initial state radiation absence) all possible sources that may lead to the 
P{' — P^ et disbalance in a final state. We use here the methods applied in [fTOl — [fT~8l (see also 
(20)) and ll2Tl . ll2~2l for analogous task at LHC energy. The corresponding cuts on physical 
variables, introduced in [fTOl — fTlTl . are applied here. Their efficiency is estimated at the particle 
level of simulation at Tevatron energy with account of DO detector geometry. 

We consider here the case of the Tevatron Run II luminosity L = 10 32 cm~ 2 s~ l . It will 
be shown below that its value is quite sufficient for selecting the event samples of large enough 
volume for application strict cuts as well as of new physical variables introduced in [ lOI- lfTTI . 

Section 2 is a short introduction into the physics connected with the discussed problem. 
General features of "7 + jet" processes are presented here. We review the possible sources of 
the P{' and P^ et disbalance and the ways of selecting those events where this disbalance has a 
minimal value on the particle level. 

In Section 3.1 we give the definitions are given for the transverse momenta of different 
physical objects that we suppose to be important for studying the physics connected with a jet 
calibration procedure. Values of these transverse momenta enter into the P r balance equation that 
reflects the total P t conservation law for the pp-collision event as a whole. 

Section 3.2 describes the criteria we have chosen to select "7 + jet" events for the jet 
energy calibration procedure. The "cluster" (or mini-jet) suppression criterion {Ptcur) which 
was formulated in an evident form in our previous publications iflOl — lfl~8l is used here 1 . (Its 
important role for selection of events with a good balance of Pt 1 and P^ et will be illustrated in 
Sections 5-8.) 2 These clusters have a physical meaning of a part of another new experimentally 
measurable quantity, introduced in lITOl - lfTHI for the first time, namely, the sum of P t of those 
particles that are out of the "7 + jet" system (denoted as p t out ) and are detectable in the whole 
pseudorapidity rj region covered by the detector 3 . The vector and scalar forms of the total P t 
balance equation, used for the pp— event as a whole, are given in Sections 3.1 and 3.3 respectively. 

Another new thing is a use of a new physical object, proposed also in ifTOl - lfTHI and named 
an "isolated jet". This jet is contained in the cone of radius R = 0.7 in the 77 — <\> space and does 
not have any noticeable P t activity in some ring around. The width of this ring is taken to be of 
AR = 0.3 (or approximately of the width of 3 calorimeter towers). In other words, we will select 
a class of events having a total P t activity inside the ring around this "isolated jet" within 3 — 5% 

'We use here, as in 1131 - 1181 . for most application the PYTHIA's default jetfinder LUCELL as well as UA1 taken 
from the CMS program of fast simulation CMS JET 1241 for defining jets in an event. 

2 The analogous third jet cut thresholds E^, (varying from 20 to 8 GeV) for improving a single jet energy resolu- 
tion in di-jet events were used in 1291 . 
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of jet P t . It will be shown in Sections 6, 7 and Appendix 2 that the number of events with such a 
clean topological structure would not be small at Tevatron energy and L = 10 32 cm~ 2 s~ 1 . 

Section 4 is devoted to the estimation of a size of a non-detectable neutrino contribution to 
a jet. The correlation of the upper cut value, imposed onto p t miss 4 5 with the mean value of P t of 
neutrinos belonging to the jet P t is considered. The detailed results of this section are presented 
in the tables of Appendix 1. They also include the ratios of the gluonic events qg — > q + 7 
containing the information about the gluon distribution inside a proton. In the same tables the 
expected number of events (at L int = 300 pb^ 1 ) having charm (c) and beauty (b) quarks in the 
initial state of the gluonic subprocess are also given. 

Since the jet energy calibration is rather a practical than an academic task, in all the follow- 
ing sections we present the rates obtained with the cuts varying from strict to weak because their 
choice would be a matter of step-by-step statistics collection during the data taking. 

Section 5 includes the results of studying the dependence of the initial state radiation (ISR) 
Pt-spectrum on the cut imposed on the clusters P t (Ptcur) an d on me an gle between the trans- 
verse momenta vectors of a jet and a photon. We also present the rates for four different types of 
"7 + J e ^" events, in which jet fits completely in one definite region of the calorimeter: in Central 
Calorimeter (CC) with \r)\ < 0.7 or in Intercryostat Calorimeter (IC) with 0.7< \r)\ < 1.8 or in End 
Calorimeter (EC) with 1.8 < \r]\ < 2.5 or, finally, in Forward Calorimeter (FC) with 2.5 < \r)\ < 4.2. 

In Section 6 our analysis is concentrated on the "7 + 1 jet" events having a jet entirely 
contained within the central calorimeter region. The dependence of spectra of different physical 
variables 5 (and among them those appearing in the P t balance equation of event as a whole) on 
Ptcur is shown there. 

The dependence of the number of events (for L int = 300 pb^ 1 ) on Ptcur as wei l as me 
dependence on it of the fractional {P{' — P t Jet )/P t y disbalance is studied in Section 7. The 
details of this study are presented in the tables of Appendix 2 that together with the corresponding 
Figs. 10-12 can serve to justify the variables and cuts introduced in Section 3. 

In Section 8 we present an estimation of the efficiency of background suppression (that was 
one of the main guidelines to establish the selection rules proposed in Section 3) for different 
numerical values of cuts. 

The importance of the simultaneous use of the above-mentioned new parameters Ptcur 
and Ptcur an( ^ a ^ so °f me "isolated jet" criterion for background suppression (as well as for 
improving the value of the Pt 1 and P t 3et balance) is demonstrated in Tables 8-1 1 of Section 8 as 
well as in the tables of Appendix 3 for various Pt' intervals. 

The tables of Appendix 3 include a fractional disbalance values (Pt 1 — P t Jet ) / Pf 1 that are 
found with an additional (as compared with tables of Appendix 2) account of the P t out cut. They 
contain the final and first main result (as they include the background contribution) of our study of 
setting an absolute scale of a jet energy at the particle level defined by generation with PYTHIA. 

Section 9 contains the second main result of our study of "7 + jet" events at Tevatron 
energy. Here we investigate a possibility of using the same sample of the topologically clean 
"7 + J e ^" events, obtained with the described cuts, for determining the gluon distribution in a 
proton (as it was done earlier for LHC energy in JT81 . Il2"0l ). The kinematic plot presented here 
shows what a region of x and Q 2 variables can be covered at Tevatron energies with a sufficient 
number of events for this aim. The comparison with the kinematic regions covered by other 

4 see (7) for definition 

5 mostly those that have a strong influence on the P t 7 — Pt 3et balance in an event. 



experiments where parton distributions were studied is also shown in the same plot (see Fig. 17). 

About the Summary. We tried to write it in a way allowing a dedicated reader, who is 
interested in result rather than in method, to pass directly to it after this sentence. 

Since the results presented here were obtained with the PYTHIA simulation, we are plan- 
ning to carry out analogous estimations with another event generator like HERWIG, for example, 
in subsequent papers. 



2. GENERALITIES OF THE "7 + jet" PROCESS. 

Useful variables are introduced here for studying their effects on the initial and final state radiation 
basing on the simulation in the framework of PYTHIA. Other effects of non-perturbative nature like pri- 
mordial parton k t effect, parton-to-jet hadronization that may lead to Pp — P^ et disbalance within the 
physical models used in PYTHIA are also discussed. 



2.1 Leading order picture. 

The idea of absolute jet energy scale setting calibration) by means of the physical process ll pp — > 
7 + jet + X" was realized many times in different experiments (see [1-8] and references therein). 
It is based on the parton picture where two partons (qq or qg), supposed to be moving in differ- 
ent colliding nucleons with zero transverse momenta (with respect to the beam line), produce a 
photon called the "direct photon". This process is described by the leading order (LO) Feynman 
diagrams shown in Fig. 1 6 for the "Compton-like" subprocess 

qg^q + i (la) 

and for the "annihilation" subprocess 

qq^g + 7- (16) 

As the initial partons were supposed to have zero transverse momenta, P t of the "7+parton" 
system produced in the final state should be also equal to zero, i.e. one can write the following P t 
balance equation for photon and final parton 

-p^y+part _ p 7 _|_ popart _ g ^2) 

One could expect that the transverse momentum of the jet produced by the final state parton (q or 
g) with p t f ar ' = — P t 7 will be close in magnitude with a reasonable precision to the transverse 
momentum of the final state photon, i.e. Pt je * ~ — Pt 7 - Thus, in principle, having a well- 
calibrated photon energy scale one can determine a jet energy scale. That is the a main idea of the 
procedure. But a more detailed analysis leads to some features needed to be taken into account 
and to a photon-jet P t balance equation in a more complex form. 




(a) (b) 
Fig. 1 : Some of the leading order Feynman diagrams for direct photon production. 
6 for the explanation of the numeration of lines see Section 2.2 



2.2 Initial state radiation. 

Since we believe in the perturbation theory, the leading order (LO) picture described above is 
expected to be dominant and determine the main contribution to the cross section. The Next-to- 
Leading Order (NLO) approximation (see some of the NLO diagrams in Figs. 2 and 4) introduces 
some deviations from a rather straightforward LO-motivated idea of a jet energy calibration. A 
gluon radiated in the initial state (ISR), as it is seen from Fig. 2, can have its own non-zero 
transverse momentum p t 9 luon = p^SR _^ q Apart of a problem of appearance of extra jets (or 
mini-jets and clusters), that will be discussed in what follows, it leads to the non-zero transverse 
momenta of partons that appear in the initial state of fundamental 2 — > 2 QCD subprocesses (la) 
and (lb). As a result of the transverse momentum conservation there arises a disbalance between 
the transverse momenta of a photon P 4 7 and of a parton p t part produced in the fundamental 2^2 
process 5 + 6^7 + 8 shown in Fig. 2 (and in Fig. 3) and thus, finally, the disbalance between 
Pt' and P t of a jet produced by this parton. 




Fig. 2: Some of Feynman diagrams of direct photon production including gluon radiation in the initial state. 



Following lfP3l - lfP7l and ll25l we choose the modulus of the vector sum of the transverse 
momentum vectors P t 5 and P t 6 of the incoming into 2^2 fundamental QCD subprocesses 
5 + 6 — > 7 + 8 partons (lines 5 and 6 in Fig. 2) and the sum of their modulus as two quantitative 
measures 

P t 5+6 = |P t 5 + P t 6 |, P t 56 = |P t 5 | + |P t 6 | (3) 
to estimate the P t disbalance caused by ISR 1 . The modulus of the vector sum 

p7+jet = |p t 7 + p t jet| (4) 

was also used as an estimator of the final state P t disbalance in the "7 + jet" system in 1T3l — ifTTl . 

The numerical notations in the Feynman diagrams (shown in Figs. 1 and 2) and in formula 
(3) are chosen to be in correspondence with those used in the PYTHIA event listing for description 
of the parton-parton subprocess displayed schematically in Fig. 3. The "ISR" block describes the 
initial state radiation process that can take place before the fundamental hard 2^2 process. 
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Fig. 3: PYTHIA "diagram" of 2 -> 2 process (5+6^7+8) following the block (3+4^5+6) of initial state radiation 
(ISR), drawn here to illustrate the PYTHIA event listing information. 

7 The variable Ft 5+6 was used in analysis in 1 1 01 — 1 1 31 . 



2.3 Final state radiation. 

Let us consider fundamental subprocesses in which there is no initial state radiation but instead 
final state radiation (FSR) takes place. These subprocesses are described in the quantum field 
theory by the NLO diagrams like those shown in Fig. |4j It is clear that appearance of an extra 
gluon leg in the final state may lead to appearance of additional jets (or clusters) in an event as it 
happens in the case of ISR described above. So, to suppress FSR (manifesting itself as some extra 
jets or clusters) the same tools as for reducing ISR should be used. But due to the string model of 
fragmentation used in PYTHIA it is much more difficult to deduce basing on the PYTHIA event 
listing information the variables (analogous to (3) and (4)) to describe the disbalance between 
P t of a jet parent parton and P t 7 . That is why, keeping in mind a close analogy of the physical 
pictures of ISR and FSR (see Figs.|2]and|4|), we shall concentrate in the following sections on the 
initial state radiation supposing it to serve in some sense as a quantum field theory perturbative 
model of the final state radiation mechanism. 




Fig. 4: Some of Feynman diagrams of direct photon production including gluon radiation in the final state. 

2.4 Primordial parton k t effect. 

Now after considering the disbalance sources connected with the perturbative corrections to the 
leading order diagrams let us mention the physical effects of a non-perturbative nature. A pos- 
sible non-zero value of the intrinsic transverse parton velocity inside a colliding proton may be 
another source of the Pi and p t part disbalance in the final state. Nowadays this effect can be 
described mainly in the phenomenological way. Its reasonable value is supposed to lead to the 
value k t < 1.0 GeV/ c. Sometimes in the literature a total effect of ISR and of the intrinsic parton 
transverse momentum is denoted by a common symbol "kt". Here we follow the approach and the 
phenomenological model used in PYTHIA where these two sources of the Pi' and Pi et disbal- 
ance, having different nature, perturbative and non-perturbative, can be switched on separately by 
different keys 8 . In what follows we shall keep the value of k t mainly to be fixed by the PYTHIA 
default value (kt) = 0.44 GeV jc. The dependence of the disbalance between Pi 1 and Pi et on 
possible variation of k t was discussed in detail in lfl7l[i"9l . The general conclusion from there is 
that any variation of k t within reasonable boundaries (as well as slightly beyond them) does not 
produce a large effect in the case when the initial state radiation is switched on. The latter makes 
a dominant contribution. 

2.5 Parton-to-jet hadronization. 

Another non-perturbative effect that leads to the Pt 1 — Pt 3et disbalance is connected with hadroniza- 
tion (or fragmentation into hadrons) of the parton produced in the fundamental 2^2 subprocess 
into a jet. The hadronization of a parton into a jet is described in PYTHIA within the Lund string 
fragmentation model. The mean values of the fractional (P t jet — pP arton y p t P arton disbalance is 
presented in the tables of Appendix 2 for UA1 jetfinding algorithm. Is is seen that a hadronization 
effect has a sizable contribution into Pi' — Pi et disbalance. 

Variables MSTP(61) for ISR and PARP(91), PARP(93), MSTP(91) for intrinsic parton transverse momentum k t 
(see 0) 



3. CHOICE OF MEASURABLE PHYSICAL VARIABLES FOR THE "7 + jet" PRO- 
CESS AND THE CUTS FOR BACKGROUND REDUCTION. 

The classification of different physical objects that participate in "7 + jet" events and that may 
give a noticeable contribution into the total Pfbalance in an event as a whole is done. 

Two new physical observables, namely, Pt of a cluster and Pt of all detectable particles beyond 
"7 + jet" system, as well as the definition of isolated jet, proposed for studying Pt 1 — Pt jet disbalance in 
IHW - M7V . are discussed. 

The selection cuts for physical observables of "7 + jet" events are given. 

The Pfbalance equation for the event as a whole is written in scalar form that allow to express the 
Pt 1 — Pt 3&t disbalance in terms of the considered physical variables. 

Apart from (la) and (lb), other QCD subprocesses with large cross sections, by orders of 
magnitude larger than the cross sections of (la) and (lb), can also lead to high P t photons and 
jets in final state. So, we face the problem of selecting signal "7 + jet" events from a large QCD 
background. Here we shall discuss a choice of physical variables that would be useful, under 
some cuts on their values, for separation of the desirable processes with direct photon ("j dir ") 
from the background events. A possible "7^" — candidate" may originate from the 7r°, 77, cu and 
K® meson decays or may be caused by a bremsstrahlung photon or by an electron (see Section 
8). 

We take the DO ECAL size to be limited by j 77 1 <2.5 and the calorimeter to be limited by 
\rj\ < 4.2 and to consist of CC, IC, EC, FC parts, where n = —In (tan (0/2)) is a pseudorapidity 
defined in terms of a polar angle 9 counted from the beam line. In a plane transverse to the beam 
line the azimuthal angle defines directions of p t Jet and P t 7 . 



3.1 Measurable physical variables and the P t vector balance equation. 

In pp — > 7 + jet + X events we are going to study the main physical object will be a high P t jet 
to be detected in the \rj\ < 4.2 region and a direct photon registered by the ECAL up to \n\ < 2.5. 
In these events there will be a set of particles mainly caused by beam remnants, i.e. by spectator 
parton fragments, that are flying mostly in the direction of a non-instrumented volume (\n\ > 4.2) 
in the detector. Let us denote the total transverse momentum of these non-observable particles (i) 
as 

£ EE Pj^- 2 . (5) 

iG|»?|M.2 

Among the particles with \rj\ < 4.2 there may also be neutrinos. We shall denote their total 
momentum as 

E Pt f M = PtM- (6) 

ie|??|<4.2 

A sum of transverse momenta of these two kinds of non-detectable particles will be denoted as 

p miss 9 . 

Pt miss = Pt M + P t l??l>4 - 2 . (7) 

A high-energy jet may also contain neutrinos that may carry a part of the total jet energy. 
The average values of this energy can be estimated from a simulation. 

9 This value is a part of true missing P t in an experiment that includes the detector effects (see [1, 2]). 



From the total jet transverse momentum P t e we shall separate the part that, in principle, 
can be detected in the ECAL+HCAL calorimeter system and in the muon system. Let us denote 
this detectable part as P t jet (small "j"!). So, we shall present the total jet transverse momentum 
P t Jet as a sum of three parts: 

1. Ptjv*, containing the contribution of neutrinos that belong to the jet, i.e. a non-detectable 
part of jet P t (i - neutrino): 

p tM=E p t(,)- (8) 

id Jet 

2. Pt(u)> containing the contribution of jet muons to P t Je< (j - muon): 

iS Jet 

These muons make a weak signal in the calorimeter but their energy can be measured, in 
principle, in the muon system (in the region of \r]\ < 2.5 in the case of DO geometry). Due to the 
absence of the muon system and the tracker beyond the \r)\ < 2.5 region, there exists a part of 
P t Jet caused by muons with \rj\ > 2.5. We denote this part as Ptj\fL| >2 5)- It is non-detectable 
part and can be considered as an analogue of Ptf^y 

As for both points 1 and 2, let us say in advance that the estimation of the average values 
of neutrino and muon contributions to P t Jet (see Section 4 and Tables 1-3 of Appendix 1) have 
shown that they are quite small: about 0.30% of (Pt Jet ) a ii is due to neutrinos and about 0.33% of 
(Pt Jet )aii is due to muons where "a//" means averaging over all events including those without 
neutrinos and/or muons in jets. So, they together may cause approximately about 0.63% of the 
PC and P^ et disbalance if muon signal is lost. 

3. Finally, as we have mentioned before, we use Pt- 76 * to denote the part of P t Jet which 
includes all detectable particles of the jet 10 , i.e. the sum of P t of jet particles that may produce a 
signal in the calorimeter (paid) and muon system (fx): 

Pt je * = Ptgio) + Pi 6 ?, |^| < 2.5. (10) 
Thus, in the general case we can write for any 77 values: 

PJet i ) jet , td Jet , -p Jet /i i \ 

t "ft + -tt(V) + ■rt( At) | J? M|>2.5)- V^J 

In the case of pp — > 7 + jet + X events the particles detected in the \r)\ < 4.2 region may 
originate from the fundamental subprocesses (la) and (lb) corresponding to LO diagrams shown 
in Fig. 1, as well as from the processes corresponding to NLO diagrams (like those in Figs. 2, 4 
that include ISR and FSR), and also from the "underlying" event [1], of course. 

So, for any event we separate the particles in the \rj\ < 4.2 region into two subsystems. 
The first one consists of the particles belonging to the "7 + Jet" system (here "Jet" denotes the 
jet with the highest P t , greater 30 GeV/c, having the total transverse momentum Y* t 1+Jet (large 
"Jet", see (4)). The second subsystem involves all other (O) particles beyond the "7 + Jet" sys- 
tem in the region, covered by the detector, i.e. \r)\ < 4.2. The total transverse momentum of this 
O-system are denoted as P t and it is a sum of P t of additional mini-jets (or clusters) and P t of 
single hadrons, photons and leptons in the \r)\ < 4.2 region. Since a part of neutrinos are also 
present among these leptons, the difference of Pt( v ) an d Ptjv* gives us the transverse momentum 

10 We shall consider the issue of charged particles contribution with small P t into the total jet P t while discussing 
the results of the full GEANT simulation (with account of the magnetic field effect) in our forthcoming papers. 



Pt&, = Pt(„) - Pt$ |i?1<4.2, (12) 

carried out by the neutrinos that do not belong to the jet but are contained in the \r) \ < 4.2 region. 

We denote by P t ou * a part of P t ° that can be measured, in principle, in the detector. Thus, 
P t out is a sum of P t of other mini-jets or, generally, clusters (with p t dust smaller than P t Jet ) and 
P t of single hadrons (h), photons (7) and electrons (e) with < 4.2 and muons (/S) with \r]^\ < 2.5 
that are out of the "7 + jet" system. For simplicity these mini-jets and clusters will be called 
"clusters" n . So, for our "7 + jet" events p t out is the following sum (all {h, 7, e, fj,} £ Jet): 

Pout _ td dust 1 -p sing -p nondir 13 ,p O . I I ^ a ry 

t -M + rt(fc) + *t( 7 ) + ^t(e) + rt(^,|^M|<2.5)) |??|<4.Z. (^IJJ 

And thus, finally, we have: 

Pt° = ~P t ° Ut + Ptfu) + Pt^i.^^2.5)- (14) 

With these notations we come to the following vector form lfi~3l of the P t - conservation law for 
the "7 + Jet" event (where 7 is a direct photon) as a whole (supposing that the jet and the photon 
are contained in the corresponding detectable regions): 

PJ + P/ et + P t ° + P t l"l >4 - 2 = (15) 
with last three terms defined correspondingly by (11), (15) and (5) respectively. 



3.2 Definition of selection cuts for physical variables. 

1. We shall select the events with one jet and one "7 dir -candidate" (in what follows we shall 
designate it as 7 and call the "photon" for brevity) 12 with 

Pf > 40 GeV/c and P t Jet > 30 GeV/c. (16) 

In the simulation the ECAL signal is considered as a candidate for a direct photon if it fits inside 
one DO calorimeter tower having size 0.1 x 0.1 in the 77 — (ft space. 

For most of our applications in Sections 4, 5 and 6 mainly the PYTHIA jetfinding algorithm 
LUCELL will be used 13 . The jet cone radius R in the 77 — space counted from the jet initiator 
cell (ic) is taken to be R ic = ((Ar?) 2 + (A0) 2 ) 1 / 2 = 0.7. 

2. To suppress the contribution of background processes, i.e. to select mostly the events with 
"isolated" direct photons and to discard the events with fake "photons" (that may originate as 
7 d * r -candidates from meson decays, for instance), we restrict 

a) the value of the scalar sum of P t of hadrons and other particles surrounding a "photon" 
within a cone of R] sol = ((Ar?) 2 + (A0) 2 ) 1 / 2 = 0.7 ("absolute isolation cut") 14 

Y. p t=Pt s ° l <PtcuT^ (17) 

i€R 



11 As was already mentioned in Introduction, these clusters are found by the LUCELL jetfinder with the same 
value of the cone radius as for jets: R clust = RJ et = 0.7. 

12 only in Section 8, devoted to the backgrounds, we shall denote 7 rfjr -candidate by 7 
13 Comparison with the UA1 and UA2 jetfinding algorithms was presented in ll9U15lfT6l 

14 We have found that S/B ratio with R] sol = 0.7 is in about 1.5 times better than with Rj sol = 0.4 what is 
accompanied by only 10% of additional loss of the number of signal events. 



b) the value of a fraction ("fractional isolation cut") 

Y^P t yP t ^^<el UT . (18) 

3. We accept only the events having no charged tracks (particles) with P t > 5 GeV/c within the 
R = 0.4 cone around the 7 Q!ir -candidate. 

4. To suppress the background events with photons resulting from tt°, r], lo and if <j meson decays, 
we require the absence of a high P t hadron in the tower containing the 7 cfer -candidate: 

P t hadr < 7 GeV/c. (19) 

At the PYTHIA level of simulation this cut may effectively take into account the imposing of 
an upper cut on the HCAL energy in the cells behind the ECAL signal cells fired by the direct 
photon. In real experimental conditions one can require that a fraction of the photon energy, 
deposited in ECAL to be greater than some value (~ 0.95 — 0.96 as it is now at DO). 

5. We select the events with the vector P t Jet being "back-to-back" to the vector P t 7 (in the plane 
transverse to the beam line) within A</> defined by the equation: 

(7 , jet) = 180° ± A0, (20) 

where (7i jet ) is the angle between the Pt 1 and P t jet vectors: P t 7 P t Jei = Pt 1 P t Jet -cos(0( 7i jet )), 
P t -y = |PJ|, P t Jet = \P t Jet \. The cases A</> < 17°, 11°, 6° are considered in this paper (6° is 
approximately one DO calorimeter tower size in 0). 

6. As we have already mentioned in Section 3.1, one can expect reasonable results of the jet 
energy calibration procedure modeling and subsequent practical realization only if one uses a set 
of selected events with small p t miss (see (7) and d25t ). So, we also use the following cut: 

r> miss ^ -r> miss /O 1 \ 

P S PcUT- \ li > 

The aim of the event selection with small p t miss is quite obvious: we need a set of events with a 
reduced P t Jet uncertainty due to a possible presence of a non-detectable particle contribution to 
a jet and due to the term p t ^> 4 - 2 (see (7) and (I25l0. 

The influence of Ptcur on tne selection of events with a reduced value of the total sum of 
neutrino contribution into P t Jet is studied in Section 4. 

7. The initial and final state radiations (ISR and FSR) manifest themselves most clearly as some 
final state mini-jets or clusters activity. To suppress it, we impose a new cut condition that was 
not formulated in an evident form in previous experiments: we choose the "7 + jet" events that 
do not have any other jet-like or cluster high P t activity by selecting the events with the values 
of P t clust (the cluster cone R c iust(Vi 0) = 0-7), being lower than some threshold P t % u v st T value, i.e. 
we select the events with 

P. dust r> dust /'00^ 
t — PcUT \ LL > 

(PtouT = 10, 5 GeV/c are most effective as will be shown in Sections 6-8). Here, in contrast 
to lfT3l — lfT7l . the clusters are found by one and the same jetfinder LUCELL while three different 
jetfmders UA1, UA2 and LUCELL are used to find the jet (P t Jet > 30 GeV/c) in the event. 

8. Now we pass to another new quantity (proposed also for the first time in lfT3l - lfT71 ) that can 
be measured at the experiment. We limit the value of the modulus of the vector sum of P t of all 
particles, except those of the "7 + jet" system, that fit into the region \r]\ < 4.2 covered by the 
ECAL and HCAL, i.e., we limit the signal in the cells "beyond the jet and photon" region by the 



following cut: 



P° ut < Ptcur, l^l<4.2. (23) 



The importance of Ptchr an d Ptcirr f° r selection of events with a good balance of PP and 
Pp et and for the background reduction will be demonstrated in Sections 7 and 8. 

Below the set of selection cuts 1-8 will be referred to as "Selection 1". The last two of 
them, 7 and 8, are new criteria [Tl~3l not used in previous experiments. 

9. In addition to them one more new object, introduced in [|T3l - ifTTll and named an "iso- 
lated jet", will be used in our analysis, i.e. we shall require the presence of a "clean enough" (in 
the sense of limited P t activity) region inside the ring of AR = 0.3 width (or approximately of a 
size of three calorimeter towers) around the jet. Following this picture, we restrict the ratio of the 
scalar sum of transverse momenta of particles belonging to this ring, i.e. 

P t Hn9 /P/ et = e jet <4 e \ where P t rmg = £ |P t *|. (24) 

ie0.7<R<1.0 

(e 3 et is chosen to be 3 — 5%, see Sections 7 and 8). 

The set of cuts 1-9 will be called in what follows "Selection 2". 

The exact values of the cut parameters PtSuT> e cuT> ejet > PtcuT> Ptcur Wlil be specified 
below, since they may be different, for instance, for various Pt 1 intervals (being looser for higher 

Pt 1 )- 



3.3 The scalar form of the P t balance equation and the jet energy calibration procedure. 

Let us rewrite the basic P t -balance equation (fl~5T) of Section 3.1 with the notations introduced here 
in the scalar form more suitable for the following applications: 

Pt 1 - Pt Jet , 



P 



7 



cosA(f)) + P t (0+r] > 4.2) /PP, (25) 



where P t (0+r] > 4.2) = (P t ° + pj^ 4 ' 2 )) . n Jet with n Jet = P t Jet /Pt • 

As will be shown in Section 7, the first term on the right-hand side of equation (l25l) . i.e. ( 1 — 
cosAcj)) is negligibly small as compared with the second term 15 and tends to decrease fast with 
growing P t Jet . So, in this case the main contribution to the P t disbalance in the "7 + jet" system 
is caused by the term P t (0 + rj> 4.2) / PP '. 

P t Jet can be easily expressed from equation (1231) through: 

P t Jet = a-PP (26) 

with a defined as a = cosAcj) - P t (0+ rj > 4.2) jPP . 

Having defined in every selected event P t Jet from equation (I26T) one can determine calibra- 
tion coefficients {Q} via minimizing of a standard deviation of the function: 




P t Jet -EZiCiPi 



AP 



Jet 



(27) 



15 in a case of Selection 1 



In this expression Ni is a number of calorimeter layers 16 , P t l,c is energy deposition in the i— th 
calorimeter layer and AP t Jet is the error on P t Jet caused by uncertainty in a (Aa) and uncertainty 
due to limited accuracy of Pt 7 determination (APt 7 ) 17 . So, one can write (see (1261)): 

AP t Jet = Aa © AP t 7 (28) 

Obtained in this way the calibration coefficients {Ci} in the selected "7 + jet" events for 
every bin of rf Gt and calorimeter P t ^ et then should be applied to energy depositions in each layer 
P/' c of a found jet in any event to reconstruct a jet transverse momentum at the particle level The 
accuracy of such a reconstruction will directly depend on the accuracy of the coefficients {Ci}. 
The latter, in their turn, are caused by the error of AP t Jet (see G3) 18 • 

Having determined relatively perfectly a photon energy scale the AP t 7 , the AP t Jet uncer- 
tainty will be mainly defined by Aa (namely by term P t (0+rj > A.2)/P{ 1 of equation (1251). 

4. ESTIMATION OF A NON-DETECTABLE PART OF P t Jet AND P t miss SPECTRA. 

It is shown that by imposing an upper cut on the missing transverse momentum p t mtss < 10 GeV/c 
one can reduce the correction to the measurable part of Pt jet due to neutrino contribution down to the 
value ofA u = (Ptffyaii events = 0.1 GeV/c in all intervals ofPf. 

In Section 3.1 we have divided the transverse momentum of a jet, i.e. P t Jet , into two parts, 
a detectable P t jet and non-detectable (P t Jet - P t jet ), consisting of P t J $ and Ptf^\ v \ >2 .B) (see (11)). 
In the same way, according to equation (15), we divided the transverse momentum P t of "other 
particles", that are out of i y dir + jet system, into a detectable part P t out and a non-detectable part 
consisting of the sum of P t f^ and PtQ i \ v \>2.5) ( see (15)) 

We shall estimate here what part of P t Jet may be carried out by non-detectable particles 
(mainly neutrinos originating from weak decays) 19 . 

We shall consider the case of switched-on decays of 7r ± and K ± mesons 20 . Here and 
K ± meson decays are allowed inside the solenoid volume with the barrel radius R B = 80 cm and 
the distance from the interaction vertex to Endcap along the z-axis L = 130 cm (DO geometry). 

For this aim we shall use the bank of the signal "7+ jet" events, i.e. caused by subprocesses 
(la) and (lb), generated for three P t 7 intervals: 40 < P t 7 < 50, 70 < P t 7 < 90 and 90 < P 7 < 
140 GeV/c and selected with conditions (16)— (23) (Selection 1) and the following cut values: 

P t ^ T = 4GeV/c, el UT = 7%, A0<17°, P t %% = 30 GeV/c. (29) 

Here the cut PtcuT = 30 GeV/c has the meaning of a very weak restriction on mini-jets or 
clusters activity. No restriction was imposed on the P t out value. The results of analysis of these 
events, based on the application of LUCELL jetfinder, are presented in Fig.|5] 

16 iV ( = 8 at DO (4 for ECAL and 4 for HCAL). 

17 For instance, in the central region of DO calorimeter (\r/\ < 0.9) electron/photon energy resolution can be written 
through a jE = 23%/y/[E) © 20%/ E © 0.4%. 

18 Other possibility, based on the usage of artificial neural networks (ANN), was also considered (see 1271 and 
1281 1. In this approach one can obtain a better energy resolution of the reconstructed jet but it requires a bigger 
statistics for ANN training. The calibration coefficients {Ci} in this case will be replaced by set of ANN weights 
{wij} and function J27t — by a more complicated expression. 

19 In Q7| and |26| it was shown that main source of high P t neutrinos in background processes are W ± decays, 
which also contain e ± that in its turn may fake direct photons. 

20 According to the PYTHIA default agreement, tt^ and mesons are stable. 



40<Pt Y <50 GeV/c 



70<Pt y < 90 GeV/c 



90<Pt Y <140 GeV/c 






Fig. 5: a) p t mias spectra in all events; b) p t mlss spectra in events having jets with non-zero P t neutrinos, i 
Pt(v) > 0; c ) Ptfu) spectra and their mean values dependence on the values of PtcuT m various P t 7 (« P t Je 
intervals. 7r ± and K ± meson decays are allowed inside the solenoid of R = 80 cm and L = 130 cm (P t 
30 GeV/c). 



The first row of Fig. |5] contains p t miss spectra in the "7 + jet" events for different Pt 1 
intervals and demonstrates (to a good accuracy) their practical independence on Pp '. 

In the second row of Fig. |5] we present the spectra of p t miss for those events (denoted as 
Pt{y) > 0) which contain jets having neutrinos, i.e. having a non-zero Ptu/j component of P t Jet . 
These figures show a very weak dependence of the P t miss spectrum on the direct photon P^ 
Comparison of the number of entries in the second row plots of Fig. |5] with those in the first row 
allows to conclude that the part of events with the jet having the non-zero neutrinos contribution 
is about 15 - 18%. 

The effect of imposing general Ptcur m eacn event of our sample is shown in the third 
row of Fig. |3 The upper cut Ptcur = 1000 GeV/c means the absence of any upper limit for 
PthA- The most important illustrative fact that in the absence of any restriction on p t m%ss the 
total neutrino P t inside the jet averaged over all events can be as large as PJ^ « 0.32 GeV/c 
at 90 < Pi 1 < 140 GeV/c GeV/c (see the right-hand plot of the third row in Fig. In the 
40 < P t < 50 GeV/c interval, we have already a very small mean value of Ptiyj equal to 
0.12 GeV/c even without imposing any Ptcur- From the same plots of the third row of Fig.|5]we 
see that with general Ptcur = 10 GeV/ c the average correction due to neutrino contribution is 
0.1 GeV /cm all three intervals of P t 7 . 

At the same time, as it was demonstrated in [17] and ll2"6ll. this cut essentially reduces the 
admixture of the e ± -events, in which e ± , mainly originating from the W ± — > e ± u weak decays, 
may fake the direct photon signal. These events are characterized by big values of p t miss (it is 
higher, on the average, by about one order of magnitude than in the signal "<y dM " + jet" events) 
that may worsen the jet calibration accuracy. 

The situation, analogous to neutrino, holds for the P t f^j contribution. 

The detailed information about the values of non-detectable PJ^ averaged over all events 
(no cut on p t m%ss was used) as well as about mean P t values of muons belonging to jets (Ptu)} 
is presented in Tables 1-3 of Appendix 1 for the sample of events with jets which are entirely 
contained in the central region of the calorimeter (\rf et \ < 0.7) and found by UA1 jetfinder. In 
these tables the ratio of the number of events with non-zero PtuA to the total number of events is 
denoted by R^nt an( l the ratio of the number of events with non-zero P t f^ to the total number 
of events is denoted by R^vent ■ 

The quantity p t miss in events with P t J $ > is denoted in these tables as P t ^ S j et and is 
given there for three P 7 intervals (40 < P t 7 < 50, 70 < Pt 1 < 90 and 90 < Pt 1 < 140) and 
P t *** = 30, 20, 15, 10, 5 GeV /c) 21 

5. EVENT RATES FOR DIFFERENT P 7 AND rf et INTERVALS. 

The number of "7 + jet" events distribution over Pj 7 and rp is shown here. It is found that in each 
interval of the APt 7 = 10 GeV jc width the rates decrease by a factor more than 2. The number of events 
with jets which transverse momentum are completely (or with 5% accuracy) contained in CC, IC, EC and 
FC regions are presented in Tables\5\^for integrated luminosity Li n t = 300 p6 _1 . 

21 Note that the values of p t mlss and Pt™ s j e t m tne Pl° ts of Fig.|3are slightly different from those of Appendix 1 
as the numbers in from Fig.[5]were found for events in the whole \rj\ < 4.2 region. 



5.1 Dependence of distribution of the number of events on the "back-to-back" angle 0( 7i jet ) 

and on P t ISR . 

The definitions of the physical variables introduced in Sections 2 and 3 allow to study a possible 
way to select the events with a good PP and P t Jet balance. Here we shall be interested to get (by 
help of PYTHIA generator and the theoretical models therein) the form of the spectrum of the 
variable P t 56 (see (3)) (which is approximately proportional to P t ISR up to the value of intrinsic 
parton transverse momentum k t inside a proton) at different values of Pf 1 . For this aim four 
samples of "7 + jet" events were generated by using PYTHIA with 2 QCD subprocesses (la) 
and (lb) being included simultaneously. In what follows we shall call these events as "signal 
events". The generations were done with the values of the PYTHIA parameter CKIN(3)(= p ™ m ) 
equal to 20, 25, 35 and 45 GeV/c in order to cover four P t 7 intervals: 40-50, 50-70, 70-90 and 
90-140 GeV/c, respectively 22 . Each sample in these Pt intervals had a size of 5 • 10 6 events. 
The cross sections for the two subprocesses were found to be as given in Table[T] 

Table 1: The cross sections (in microbarns) of the qg — > q + 7 and qq — > g + j subprocesses for four P t 7 intervals. 



Subprocess 
type 


p r n (GeV/c) 


20 


25 


35 


45 


qg^q + 1 


0.97-10" 2 


4.78-10~ 3 


1.36-10- 3 


4.95-10" 4 


QQ^ 9 + 1 


0.20- 10" 2 


0.96-10" 3 


0.35-10" 3 


1.56-10" 4 


Total 


1.17-10" 2 


5.75-10" 3 


1.7M0" 3 


6.51-10" 4 



For our analysis we used Selection 1 (see (16)— (23)) and the values of cut parameters (32). 

In Fig.|6]we present the P t 56 spectra for two most illustrative cases of Pt 1 intervals 40 < 
Pi 1 < 50 GeV/c (two upper plots) and 70 < Pt' < 90 GeV/c (two bottom plots). The distributions 
of the number of events for the integrated luminosity L int = 300 pb^ 1 in different P t 56 intervals 
and for different "back-to-back" angle intervals 0( 7i j et \ = 180° ± A0 (A0 < 17° and 6° as well 
as without any restriction on A0, i.e. for the whole <p interval A0 < 180°) 23 are given there. 
The LUCELL jetfinder was used for determination of jets and clusters 24 . Left column of Fig. [S] 
correspond to P t dust < 30 GeV/c and serve as an illustration since it is rather a weak cut condition, 
while right column of Fig. |6] correspond to a more restrictive selection cut Ptcur = 5 GeV/c. 

Tables El and |3 show the number of events (at L int = 300 pb~ l ) left after application dif- 
ferent cuts on the angle A0 for two values of Ptcur- m me case °f weak restriction p t dust < 
30 GeV/c we can see that for the 40 < Pt 1 < 50 GeV/c interval about 75% of events are 
concentrated in the A0 < 17° range, while 41% of events are in the A0 < 6° range. As for 
70 < Pt 1 < 90 GeV/c: about 86% of events have A0 < 17° and 50% of them have A0 < 6°. 

A tendency of distributions of the number of signal "7 + jet" events to be very rapidly 
concentrated in a rather narrow back-to-back angle interval A0 < 17° as P^ grows becomes 
more distinct with a more restrictive cut on the cluster P t . From Table we see that in the first 
interval 40 < P{' < 50 GeV/c more than 99% of the events, selected with P t cur = 5 GeVjc, 
have A0 < 17°, while 76% of them are in the A0 < 6° range. It should be mentioned that after 
application of cut Ptcur = 5 GeV/c only about 40% of events remain as compared with a case 
of Pt C uT = 30 GeV/c. For 70 < P t 7 < 90 GeV/c more than 90% of the events, subject to the 
cut PtcuT = 5 GeV/c, have A0< 6°. It means that while suppressing cluster or mini-jet activity 

22 {k t ) was taken to be fixed at the PYTHIA default value, i.e. (k t ) = 0.44 GeV/c 
23 The value A<p = 6° approximately coincides with one DO HCAL tower size in the (/>-plane. 
24 More details connected with UA1 jetfinder application can be found in Section 7 and Appendix 2 for a jet 
contained in CC region. 




Fig. 6: A dependence of the number of events (at Li nt = 300 pb v ) on A(j> ma x and Ptcur f° r two Pt 1 intervals: 
40 < Ft 1 < 50 GeV/c for two upper plots and 70 < P t 7 < 90 GeV/c for two bottom plots. 



Table 2: A dependence of the number of events on A<j> max and on P t 7 for Li nt = 300 pb 

P t ^l% = 30 GeV/c. 



P<56 
(GeV/c) 




180° 


17° 


11° 


6° 


40-50 


110691 


82913 


68921 


44830 


50-70 


71075 


55132 


45716 


29692 


70-90 


14853 


12727 


10919 


7418 


90 - 140 


5887 


5534 


4974 


3655 



Table 3: A dependence of the number of events on A<j> max and on Pt 1 for Li„ t = 300 pb 

Pt^UT = 5 GeV/c. 



P t 56 
(GeV/c) 


A<fi m ax 


180° 


17° 


11° 


6° 


40-50 


37576 


37235 


35473 


27025 


50-70 


19056 


19017 


18651 


15149 


70-90 


3773 


3773 


3755 


3387 


90 - 140 


1525 


1525 


1524 


1468 



by imposing Ptcirr = 5 GeV/c we can select the sample of events with a clean "back-to-back" 
(within 17° in 0) topology of 7 and jet orientation 25 . 

Thus one can conclude that PYTHIA simulation predicts that at Tevatron energies most of 
the "7 + jet" events (more than 75%) may have the vectors P t 7 and Pt^ 6 * being back-to-back 
within A0 < 17° after imposing P t %% = 30 GeV/c. The cut P t %% = 5 GeV/c significantly 
improves this tendency 26 . 

It is worth mentioning that this picture reflects the predictions of one of the generators based 
on the approximate LO values for the cross section. It may change if the next-to-leading order or 
soft physics 27 effects are included. 

From Fig. |6] one can see that in the case when there are no restrictions on p t clust the P t 56 
spectrum becomes a bit wider for larger values of PP '. 

At the same time, one can conclude from the comparison of left and right upper plots that 
the width of the most populated part of the P t 5Q (or P/ ) spectrum reduces by about 40% with 
restricting Piqut- So, for A(p max = 17° we see that it drops from 0<P t 56<20 GeV/c for 
P t $fft = 30 GeV/c to a narrower interval of 0<P56<10 - 12 GeV/c for the = 5 GeV/c. 

At higher P t 7 intervals (two bottom plots) for the same value A(p max = 17° the reduction factor 
of the P t 56 spectrum width is more than two (from < P t 56 < 30 GeV/c for P t %% = 30 GeV/c 
to < P56 < 12 - 15 GeV/c for P£l% = 5 GeV/c). 



Thus, we can summarize that the PYTHIA generator predicts an increase in the P t 



ISR 



spectrum with growing P t 7 and this increase can be substantially reduced by imposing a restrictive 
cut on p t dust 28 . Unfortunately, P t ISR cannot be completely suppressed by A</> and P t dust cuts 
alone 29 . That is why we prefer to use the P t balance equation for the event as a whole (see 
equations (fT3l and (l23t of Sections 3.1 and 3.3), i.e. an equation that takes into account the 
ISR and FSR effects, rather than balance equation (2) for fundamental processes (la) and (lb) as 
discussed in Section 2.1 30 . 



5.2 P{' and rp dependence of event rates. 

Here we shall present the number of events for different P t 7 and 77 7 intervals as predicted 
by PYTHIA simulation with weak cuts defined mostly by (32) with only change of PtcjjT vame 
from 30 to 10 GeV/c. The lines of Table |U correspond to P 4 7 intervals and the columns to rf 1 
intervals. The last column of this table contains the total number of events (at L int = 300 pb~ l ) in 
the whole ECAL 7/ 7 -region |?7 7 | < 2.5 for a given P t 7 interval. We see that the number of events 
decreases fast with growing P t 7 (by more than 50% for each subsequent interval). 



5.3 Estimation of "7 + jet" event rates for different calorimeter regions. 

Since a jet is a wide-spread object, the r^ 6 ' dependence of rates for different P{ ( intervals will 
be presented in a different way than in Section 5.2. Namely, Tables |5]-|51 include the rates of 

25 Unfortunately, as it will be discussed below and is seen in Fig.|6] it does not mean that PtouT a U° ws to suppress 
completely the ISR. (see also the event spectra over p t clust in Fig. 7 of the following Section 6.) 

26 A growth of P t 7 produces the same effect, as is seen from Tables|2]and|3]and will be demonstrated in more detail 
in Section 6 and Appendix 2. 

27 We thank E. Pilon and J. Ph. Jouliet for the information about new Tevatron data on this subject and for clarifying 
the importance of NLO corrections and soft physics effects. 

28 for more details see Sections 6 and 7 

29 In Sections 7, 8 the effect of the additional PtcvT w ^ ^ e discussed. 

30 In Section 6 we shall study a behavior of each term that enter equation (125 \ in order to find the criteria that would 
allow to select events with a good balance of P t 7 and Pt jet . 



Table 4: Rates for L mt = 300 pb- 1 for different P t 7 intervals and rf (Ptcur = 10 GeV/c and < 17°). 



P 7 

(GeV/c) 


|?7 7 intervals 


all 


0.0-0.4 


0.4-0.7 


0.7-1.1 


1.1-1.4 


1.4-1.8 


1.8-2.1 


2.1-2.5 


0.0-2.5 


40-50 


10978 


11232 


10604 


10337 


9662 


8051 


5806 


66679 


50-60 


4483 


4210 


4489 


3938 


3624 


2814 


1562 


25121 


60-70 


2028 


1732 


1890 


1587 


1442 


984 


607 


10270 


70-80 


949 


931 


937 


753 


637 


392 


170 


4770 


80-90 


508 


513 


469 


363 


309 


180 


62 


2405 


90 -100 


302 


287 


252 


201 


149 


80 


25 


1295 


100-120 


285 


280 


257 


189 


125 


61 


11 


1207 


120-140 


134 


121 


98 


63 


38 


9 


1 


465 


40 -140 


19662 


19302 


18992 


17427 


15986 


12571 


8245 


112216 



Table 5: Selection 1. AP t jet / P t 3et = 0.00 (Ptcur = 10 GeV/c, < 17° and L mt = 300 pb' 1 ). 





cc 


CC^IC 


IC 


IC^CC, EC 


EC 


EC^IC, FC 


FC 


FC^EC 


40-50 


9965 


13719 


8152 


22225 


617 


8854 


554 


1912 


50-60 


4009 


5597 


3104 


8791 


207 


2766 


109 


413 


60-70 


1754 


2515 


1339 


3615 


71 


979 


14 


93 


70-80 


930 


1195 


651 


1593 


21 


348 


1 


23 


80-90 


503 


596 


328 


811 


9 


136 





6 


90 - 100 


283 


352 


165 


421 


3 


59 





1 


100 - 120 


263 


351 


137 


389 


2 


37 








120 - 140 


118 


143 


50 


142 


1 


7 








40 - 140 


17822 


24462 


13927 


37988 


930 


13184 


678 


2448 



Table 6: Selection 1. AP t jet / P t 3et < 0.05 (Ptcur = 10 GeV/c, < 17° and L mt = 300 pb- 1 ). 





CC 


CC^IC 


IC 


IC^CC, EC 


EC 


EC^IC, FC 


FC 


FC^EC 


40-50 


17951 


5733 


20631 


9746 


4174 


5296 


1280 


1186 


50-60 


7466 


2141 


8313 


3583 


1403 


1570 


253 


269 


60-70 


3405 


863 


3553 


1401 


492 


558 


39 


68 


70-80 


1699 


426 


1667 


577 


179 


190 


6 


17 


80-90 


902 


197 


838 


301 


75 


71 


3 


4 


90 -100 


528 


107 


440 


146 


31 


31 








100 -120 


537 


98 


384 


142 


19 


20 








120 -140 


223 


37 


143 


48 


5 


3 








40 -140 


32701 


9603 


35971 


15943 


6377 


7738 


1582 


1545 



events (L int = 300 pb x ) for different rf et intervals, covered by the central, intercryostat, end and 
forward (CC, IC, EC and FC) parts of the calorimeter and for different P* 7 (~ Pt Jet ) intervals. 

No restrictions on other parameters are used. The first column of Table 5 CC gives the 
number of events with the jets (found by the LUCELL jetfinding algorithm of PYTHIA), all 
particles of which are comprised entirely (100%) 31 in the CC part and there is a 0% sharing of 
Pj et (AP t jef = 0) between the CC and the neighboring IC part of the calorimeter. The second 
columns of the tables CC — > IC contain the number of events in which P t of a jet is shared 
between the CC and IC regions. The same sequence of restriction conditions takes place in the 
next columns. Thus, the IC, EC and FC columns include the number of events with jets entirely 
contained in these regions, while the EC — > IC, FC column gives the number of events where the 
jet covers both the EC and IC or EC and FC regions. From these tables we can see what number of 
events can be, in principle, most suitable for the precise jet energy absolute scale setting, carried 
out separately for the CC, EC and FC parts of the calorimeter in different PP intervals. 

The selection cuts are as in (32) but Ptcur = 10 GeV/c 

Less restrictive conditions, when up to 5% of the jet P t are allowed to be shared between 
the CC, EC and FC parts of the calorimeter, is given in Tabled Tables[5]and[6]correspond to the 
case of Selection 1 32 . 

From last summarizing line of Table |5] we see that for the entire interval 40 < PC < 
140 GeV/c PYTHIA predicts around 18000 events for CC and around 1000 events for EC at 
integrated luminosity L int = 300 pb -1 . 

6. INFLUENCE OF THE P t %% PARAMETER ON THE PHOTON AND JET P t BAL- 
ANCE AND ON THE INITIAL STATE RADIATION SUPPRESSION. 

The influence of PtcuT Parameter (defining the upper limit on Pt of clusters or mini-jets in the 
event) on the variables characterizing the Pp — P^ et balance is studied. 

In this section we shall study the specific sample of events considered in the previous section 
that may be most suitable for the jet energy calibration in the CC region, with jets entirely (100%) 
contained in this region, i.e. having 0% sharing of P^ et 33 with IC. Below we shall call them "CC- 
events". The PC spectrum for this particular set of events for Ptcur = 10 GeV/c was presented 
in the first column (CC) of Table Here we shall use three different jetfinders, LUCELL from 
PYTHIA and UA1, UA2 from CMS JET E3|. The P t dust distributions for generated events found 
by the all three jetfinders in two Pf intervals, 40 < Pi' < 50 GeV/c and 70 < PP < 90 GeV/c, 
are shown in Fig. 7 for P t ^% = 30 GeV/c and A0 < 17°. 

It is interesting to note an evident similarity of the p t dust spectra with P t 5Q spectra shown 
in Fig. 6 (see also Figs. 8, 9), what support our intuitive picture of ISR and cluster connection 
described in Section 2.2. 

Here we shall study in more detail correlation of p t clust with p t ISR mentioned above. The 
averaged value of intrinsic parton transverse momentum will be fixed at (k t ) = 0.44 GeV/c. 

The banks of 1-jet "7+ jet" events gained from the results of PYTHIA generation of 5- 10 6 
signal "7 + jet" events in each of four P t 7 intervals (40-50, 50-70, 70-90, 90-140 GeV/c) 34 

31 at the particle level of simulation! 

32 The cost of passing to Selection 2 (defined in Section 3.2 with e^ et < 3%) is a reduction of the number of events 
by factor equal to 2. 

33 at the PYTHIA particle level of simulation 
34 they were discussed in Section 5 



will be used here. The observables defined in Sections 3.1 and 3.2 will be restricted here by cuts 
of Selection 1 (16)-(23) and the cut parameters defined by (32). 

We have chosen two of these intervals to illustrate the influence of the Ptcur parameter on 
the distributions of physical variables, that enter the balance equation (l25l) . These distributions 
are shown in Fig. [5] (40 < P t 7 < 50 GeV/c) and Fig. (70 < Pf < 90 GeV/c). In these 
figures, in addition to three variables P t 56, P f l^ 4 - 2 , P t ° ut , already explained in Sections 2.2, 3.1 
and 3.2, we present distributions of two other variables, P t (0+r) > 4.2) and (1 — cosA0), which 
define the right-hand side of the P t 7 — P t jet balance equation d25t . The distribution of the 7-jet 
back-to-back angle A0 (see (22)) is also presented in Figs.[8l|9j 




The ISR describing variable P t 56 (defined by (3)) and both components of (see (I231)). (1 — 
cos A(f>) and P t (O+q > 4.2) / P t 7 , as well as P t out and A0, show a tendency to become smaller (the 
mean values and the widths) with the restriction of the upper limit on the P t dust value (see Figs.[5J 
0). It means that a jet energy calibration accuracy may increase with decreasing Ptcirr* which 
justifies the intuitive choice of this new variable in Section 3. The origin of this improvement 
becomes clear from the P t 56 density plot, which demonstrates a decrease of P t 56 (or P t ) 
values with decrease of Ptcur- m Section 2.3 we gave arguments why it may also influence FSR. 

Comparison of Fig. [8] (for 40 < P 7 < 50 GeV/c) and Fig. (for 70 < P 7 < 90 GeV/c) 
also shows that the values of A0 as a degree of back-to-backness of the photon and jet P t vectors 
in the 0-plane decreases with increasing P t 7 . At the same time P° ut and p t ISR distributions 
become slightly wider. It is also seen that the P t M >4 - 2 distribution practically does not depend on 

P 7 and P t dust 35 . 

It should be mentioned that the results presented in Figs. [8] and |9] were obtained with the 
LUCELL jetfinder of PYTHIA 36 . 



see also Appendix 2 and Fig. 2 
36 The results obtained with all jetfinders and P t 7 — Pt jet balance will be discussed in Section 7 in more detail. 






Fig. 8: LUCELL algorithm, A0 < 17°; 40 < P t 7 < 50 GeV/c. Selection 1. 






Fig. 9: LUCELL algorithm, A0 < 17°; 70 < P t 7 < 90 GeV/c. Selection 1. 



7. DEPENDENCE OF THE P-DISBALANCE IN THE "7 + jet" SYSTEM ON P t c ^j% 
and PtcuT PARAMETERS. 

It is shown that with Selection 2 ( that leads to about twice reduction of the number of events N event 
for < 70 GeV/c and to about 30 — 40% loss of them at P t 7 > 70 GeV/c) one can select at the particle 
level the events with a value of the fractional (Pp — P t Jet ) / Pp disbalance better than 1%. The number 
of events (at Li n t = 300 pb" 1 ) and other characteristics of "7 + jet" events are presented in tables of 
Appendix 2 for interval 40 < P t 7 < 140 GeV/c. 

Earlier we have introduced physical variables for studying "7 + jet" events (Section 3) and 
discussed what cuts for them may lead to a decrease in the disbalance of P t 7 and P t Jet . One can 
make these cuts to be tighter if more events would be collected during data taking. 

Here we shall study in detail the dependence of the P t disbalance in the "7 + jet" system 
on PtcuT and P°cut va l ues - For this aim we shall use the same samples of events as in Section 
5 that were generated by using PYTHIA with 2 QCD subprocesses (la) and (lb) and collected to 
cover three P t 7 intervals: 40-50, 70-90, 90-140 GeV/c. 

The normalized event distributions over (P t 7 — P t Jet )/ P 4 7 for two most illustrative P 4 7 in- 
tervals 40 < P t 7 < 50 and 70 < P 7 < 90 GeV/c are shown for a case of A0 < 17° in Fig. EH in 
different plots for three jetfinders LUCELL, UA1 and UA2. These plots demonstrate the depen- 
dence of the mean and mean square deviations on PtcuT value. 

More details on Ptcur dependence of different important features of "7 + jet" events are 
presented in tables of Appendix 2. They include the information about a topology of events and 
mean values of most important variables that characterize P t J — P t Jet disbalance. 

This information can be useful as a model guideline while performing jet energy calibra- 
tion procedure and also may serve for fine tuning of PYTHIA parameters while comparing its 
predictions with the collected real data. 

Appendix 2 contains the tables for events with P t 7 varying within three intervals: 40 < 
P 7 < 50, 70 < P 7 < 90 and 90 < P 7 < 140 GeV/c. A0 is limited there by A0 < 17°. Tables 
1-3 correspond to the events passed Selection 1 with a jet found by UA1 algorithm. Tables 4-6 
correspond to the events passed Selection 2. The latter allows to select events with the "isolated 
jet", i.e. events with the total P t activity in the AR = 0.3 ring around the jet not exceeding 3% of 
jet P t (see Section 3.2) 37 . 

The columns in all tables correspond to five different values of cut parameter Ptcur = 
30, 20, 15, 10 and 5 GeV/c. The upper lines contain the expected numbers N event of "CC 
events" (i.e. the number of signal "7 + jet" events in which the jet is entirely fitted into the CC 
region of the calorimeter; see Section 5) for the integrated luminosity L int = 300 

In the next four lines of the tables we put the values of Pt56, A0, P° ut and pj 77 ^ 4 - 2 defined 
by formulae (3), (22), (24) and (5) respectively and averaged over the events selected with a 
chosen P t cirr value. 

From the tables we see that the values of P t 56, Acj), P° ut decrease fast with decreasing 
PtcuT > while the averaged values of pj 71 ^ 4 - 2 show very weak dependence on it (practically con- 
stant). 

The following three lines present the average values of the variables (P t 7 — P t part )/Pt', 
(P t J -P t part )/P t J , (P i 7 -P/)/P i 7 (here J=Jet) that serve as measures of the P t disbalance in 

37 In contrast to the case of LHC energies, where we required in Selection 2 e jet < 6 — 8% for 40 < P t 7 < 50 
(see [25 1), at Tevatron energies, due to less P t activity in the space beyond the jet, one can impose the tighter cut 

e jet < 3%. 



the "7 + parton" and "7 + jet" systems as well as a measure of the parton-to-hadrons (Jet) 
fragmentation effect. 

The lines 9, 10 include the averaged values of P t (0 + 1] > 4.2) jPt and (1 — cos(A0)) 
quantities that appear on the right-hand side of P t 7 — P^ et balance equation d25t . 

After application of cut Acj) < 17° the value of (1 — cos(A0)) becomes smaller than the 
value of (P t (0 + r) > 4.2)/P t 7 ) in the case of Selection 1 and tends to decrease faster with grow- 
ing energy. So, we can again conclude that the main contribution into the P t disbalance in the 
"7 + jet" system, as defined by equation (l25l) . comes from the term P t (0+7] > 4.2)/ 'P t ~" '. With 
Selection 2 the contribution of (P t (0+rj > 4.2) / Pf) reduces with growing p t dust to the level of 
(1 — cos(A(p)) and even to smaller values. 

We have estimated separately the contributions of two terms P t ° • n Jet and p t M >4 - 2 . n Jet 
that enter P t (0+n > 4.2) (see (I25K 

Firstly from tables it is easily seen that P t ™ >4: - 2 has practically the same value in all P t J 
intervals and it does not depend neither on A0 nor on p t dust values being equal to 2 GeV/c up to 
a good precision 38 . 

A mean value of |p t N >4 - 2 . n Jet | contribution does not exceed ~ 0.15 GeV/c and a width 
(RMS) of the corresponding distribution contributes only 11 — 12% to the total width of the 
P t (0 + r) > 4.2) distribution. So, a mean and a width of P t (0 + r) > 4.2) are caused mainly by 
measurable term P t ■ n Jet 39 . Below in this section the cuts on the value of P t out is applied to 
select events with better Pt 1 and P^ et balance. 

The following two lines contain the averaged values of the standard deviations a(Db[y, J]) 
and a(Db[j,part\) of (P 7 - P t J )/P^{= Db[y, J]) and (P 7 - P t paH )/PP(= Db[^,part}) re- 
spectively. These two variables drop approximately by about 50% (and even more for Pi 1 > 
70 GeV/c) with restricting from P t cl %% = 30 GeV/c to 5 GeV/c for all P 7 intervals. 

The last lines of the tables present the number of generated events left after cuts. 

Two features are clearly seen from these tables 40 : 

(1) in events with A0 < 17° the fractional disbalance on the parton-photon level (Pt 1 — P t part ) / P t 7 
reduces to about 1% (or even less) after imposing P t dust < 10 GeV/c. It means that PtcuT = 
10 GeV/c is really effective for ISR suppression as it was supposed in Section 3.1. 

(2) parton-to-jet hadronization/fragmentation effect, that includes also FSR, can be estimated by 
the value of the following ratio (P t J — P t part ) / P t J . It always has a negative value. It means 
that a jet loses some part of the parent parton transverse momentum p t part . it is seen that in 
the case of Selection 1 this effect gives a big contribution into P 4 7 and P^ et disbalance even 
after application of P t %$, = 10 GeV/c. The value of the fractional {P t J -P t part ) / P t J disbal- 
ance does not vary strongly with Ptcur m me cases of Selections 2 and 3. 

We also see from the tables that more restrictive cuts on the observable p t clust lead to a 
decrease in the values of P 4 56 variable (non-observable one) that serves, according to (3), as a 
measure of the initial state radiation transverse momentum Pt ISR , i.e. of the main source of the 
P t disbalance in the fundamental 2^2 subprocesses (la) and (lb). Thus, variation of PtcuT 
from 30 GeV/c to 5 GeV/c (for A0 < 17°) leads to suppression of the P t 56 value (or P t ISR ) 
approximately by 40% for 40 < P 7 < 50 GeV/c and by w 60% for P 7 > 90 GeV/c. 

38 Let us emphasize that it is a prediction of PYTHIA. 

39 A contribution of Pw^and Pt 1^ | >2 5) ( see tE) ) m the selected event samples is a negligibly small. 
40 As was shown in I19II15I a transition from A</> < 180° to A<fi < 17° supposed to be most effective in low P t 7 
intervals, does not affect the (Pf — Pt Jet ) / 'Pt 1 disbalance strongly as compared with "jet isolation" criterion or cut 

on P t clust and P t out 



In the first two intervals with P t 7 < 90 GeV j c the decrease in Ptcifr leads to some decrease 
in the (P 7 -P J )/P 7 ratio (see Tables 1,2 of Appendix 2 and Fig. 10). In the case of 90 < P 7 < 
140 GeV/c the mean value of (P t 7 -P/)/P t 7 drops from 4.2% to 1.1% (see Table 3 of Appendix 
2). After we pass to Selection 2 (Table 4-6 of Appendix 2) this disbalance becomes of the 1% 
level and smaller but at the cost of statistics loss (by about 40 — 60%). Tables 4-6 clearly show 
the prediction of PYTHIA about the best level of jet calibration precision that can be achieved 
after application of Selection 2. 

Thus, to summarize the results presented in tables of Appendix 2, we want to underline that 
only after imposing the jet isolation requirement (see Tables 4-6 of Appendix 2) the mean values 
of Pp and P t Jet disbalance, i.e. (P t 7 — P i J )/P t 7 , for all Pf 1 intervals are contained inside the 
1% window for any p t clust < 20 GeV/ a The reduction of P t clust leads to lower values of mean 
square deviations of the photon-parton Db[y,part} and of photon-jet Db[y, J] balances. 

The Selection 2 (with Ptci/T = 10 GeV jc, for instance) leaves after its application the 
following number of events with jets entirely contained (see Section 5) in the CC region at L int = 
300 pft- 1 : 

(1) about 4000 for 40 < P 7 < 50 GeV/c, (2) about 3000 for 50 < P 7 < 70 GeV/c, 
(3) about 850 for 70 < P 7 < 90 GeV/c and (4) about 500 for the 90 < P 7 < 140 GeV/c. 

So, we can say that Selections 2, besides improving the P t 7 — P/ e * balance value, is also 
important for selecting events with a clean jet topology and for rising the confidence level of a jet 
determination. 

Up to now we have been studying the influence of the Ptcur parameter on the balance. Let 
us see, in analogy with Fig. [T0l what effect is produced by PtcuT variation 41 . 

If we P t °^ T = 5 GeV/c, keeping P t dust practically unbound by P t % u v % = 30 GeV jc, then, 
as can be seen from Fig.HU the mean and RMS values of the (P 7 — P t J )/P t 7 in the case of the 
LUCELL algorithm for 40 < P t 7 < 50 GeV/c decrease from 3.6% to 1.3% and from 14.5% to 
7.1%, respectively. For 70 < P t 7 < 90 GeV jc the mean and RMS values drop from 4.5% to 
0.7% and from 11.5% to 3.7% respectively. From these plots we also may conclude that variation 
of 

P°cut improves the P t -disbalance, in fact, almost in the same way as the variation of PtciiT- 
It is not surprising as the cluster P t activity is a part of the P t out activity. 

The influence of the Ptcur variation on the distribution of (Pt 1 — P t J )/Pt 1 is shown in 
Fig. H21 for Selection 1 with the fixed value Ptcur = 10 GeV j c. In this case the mean value 
of (P t 7 -P J )/P 7 drops from 3.2% to 1.3% for LUCELL and from 2.7% to 1.3% for UA2 
algorithms for the 40 < Pf 7 < 50 GeV j c interval. At the same time RMS value changes from 
12% to 7% for all algorithms. For interval 70 < P t 7 < 90 GeV jc the mean value of fractional 
disbalance (P 7 -P/) /P 7 decrease to to 1.2 - 1.4% at P tc % T = 10 GeV/c and to less then 1% 
at PtcuT = ^ GeV jc. Simultaneously, RMS decreases to about 3.7% for all three jetfinders. 

More detailed study of Ptcur influence on the (Pj 7 — P/ ei ) / P 7 disbalance will be continued 
in the following Section 8 (see also Appendix 3). 

So, we conclude basing on the analysis of PYTHIA simulation (as a model) that the new cuts 
Ptcur an d Ptcur introduced in Section 3 as well as introduction of a new object, the "isolated 
jet", are found as those that may be very efficient tools to improve the jet calibration accuracy. 
Their combined usage for this aim and for the background suppression will be a subject of a 
further more detailed study in Section 8. 



41 This variable enters into the expression P t (0 + rj > Aty/Pf, which makes a dominant contribution to the 
right-hand side of P t balance equation (??), as we mentioned above. 



40<Pt Y <50 GeV/c 



70<Pf<90 GeV/c 





Fig. 10: A dependence of (Pf - P t J )/Pt y on PtouT for LUCELL, UA1 and UA2 jetfinding algorithms and two 
intervals of P t 7 . The mean and RMS of the distributions are displayed on the plots. A</> < 17°. p° ut is not limited. 
Selection 1. 



40<Pt y <50 GeV/c 



70<Pf<90 GeV/c 






Fig. 11: A dependence of (Pf - P t J )/Pt y on P t °cuT for LUCELL, UA1 and UA2 jetfinding algorithms and two 
intervals of P t 7 . The mean and RMS of the distributions are displayed on the plots. A0 < 17°, Ptcur = 
30 GeV/c. Selection 1. 



40<Pt y <50 GeV/c 



70<Pf<90 GeV/c 






Fig. 12: A dependence of (Pf - P t J )/Pt y on P t °cuT for LUCELL, UA1 and UA2 jetfinding algorithms and two 
intervals of P t 7 . The mean and RMS of the distributions are displayed on the plots. A0 < 17°, PtauT = 
10 GeV/c. Selection 1. 



8. ESTIMATION OF BACKGROUND SUPPRESSION CUTS EFFICIENCY. 



The relative efficiency of "hadronic " cuts that are added to "photonic " ones, used to suppress the 
background in the case of inclusive photon measurement, is estimated at the particle level. It is shown 
that an imposing of new cuts on Pt of "clusters" (PtcfjT) an d on Pt activity in the region out of "7 + 
jet" system (PtcuT^ an d a ^ so an application of "jet isolation" criterion would allow to achieve further 
(after "photonic" cuts) fourteen-fold background suppression at the cost of four-fold loss of the signal 
"rydir _|_ j e £» even f S 

It is also shown that the imposing of PtcuT> PtcUT together with a usage of jet isolation criterion 
would lead to a substantial improvement of Pp — Pt 3et balance. 

The potentially dangerous role of a new source of background to the signal "j dir + jet" events 
caused by hard bremsstrahlung photons ("7 — brem") is demonstrated. It is shown that at Tevatron 
energy this new irreducible background may be compatible at low Pi 1 intervals with the 7r° contribution 
and it may grow faster with Pp increasing than the latter one. 



To estimate an efficiency of the selection criteria proposed in Section 3.2 we carried out the 
simulation 42 with a mixture of all QCD and SM subprocesses with large cross sections existing 
in PYTHIA 43 . The events caused by this set of the subprocesses may give a large background to 
the "7 dir + jet" signal events defined by the subprocesses (la) and (lb) 44 that were also included 
in this simulation. 

Three generations with the above-mentioned set of subprocesses were performed. Each of 
them was done with a different value of pj am parameter 45 that defines a minimal value of P t 
appearing in the final state of a hard 2^2 parton level fundamental subprocess in the case of 
ISR absence. These values were p™ m = 40,70 and 100 GeV/c. By 40 million events were 
generated for each of p™ m value. The cross sections of the above-mentioned subprocesses define 
the rates of corresponding physical events and thus appear in simulation as weight factors. 

We selected "7 dir -candidate +1 jet" events containing one 7* r -candidate (denoted in what 
follows as 7) and one jet, found by LUCELL, with P t 3et > 30 GeV/ c. Here and below, as we 
work at the PYTHIA particle level of simulation, speaking about the r ) dtr -candidate we actually 
mean, apart from r ) dir , a set of particles like electrons, bremsstrahlung photons and also photons 
from neutral meson decays that may be registered in one DO calorimeter tower of the A77 x A0 = 
0.1 x 0.1 size. 

Here we consider a set of 17 cuts that are separated into 2 subsets: 6 "photonic" cuts and 1 1 
"hadronic" ones. The first subset consists of the cuts used to select an isolated photon candidate 
in some Pt' interval. The second one includes the cuts connected mostly with jets and clusters 
and are used to select events having one "isolated jet" and limited P t activity out of "7 + jet" 
system. 

The used cuts are listed in Table |7] To give an idea about their physical meaning and 
importance we have done an estimation of their possible influence on the signal-to-background 
ratios S/B. The letter were calculated after application of each cut. Their values are presented in 
Table [8] for a case of the most illustrative intermediate interval of event generation with p™ m = 
70 GeV I c. In this table the number in each line corresponds to the number of the cut in Table 
|7] Three important lines of Table |S] are darkened because they will be often referenced to while 

42 PYTHIA 5.7 version with default CTEQ2L parameterization of structure functions is used here. 
43 ISUB=1, 2, 1 1-20, 28-31, 53, 68 (in notations of PYTHIA) 

44 ISUB=29 and 14 in PYTHIA. A contribution of another possible NLO channel gg -> 57 (ISUB=1 15) was found 
to be still negligible even at Tevatron energies. 
45 CKIN(3) in PYTHIA 



discussing the following Tables I914TT1 

The efficiencies Effsm) (with their errors) in Table[8]are defined as a ratio of the number 
of signal (background) events that passed under a cut (1-17) to the number of the preselected 
events (1st cut of this table). 



Table 7: List of the applied cuts (will be used also in Tablesl8l 41 111 . 



1. 


a) Pt 1 > 40 GeV/c, b) P t jet 


> 30 GeV/c 


9. A(p < 17°; 




c) \rf\ < 2.5, d) P t hadr 


<7GeV/c* 


; 10. p t miss /pf'< 0.10; 


2. 


P t is ° l < 5 GeV/c, < 15%; 




11. P t dust < 20 GeV/c; 


3. 






12. P t clust < 15 GeV/c; 


4. 


Pt% 1 <lGeV/c**; 
P™° r < 2 GeV/c, < 5%; 




13. P t clust < 10 GeV/c; 


5. 




14. P t out < 20 GeV/c; 


6. 


Njet < 3; 




15. P t out < 15 GeV/c; 


7. 


Njet < 2; 




16. P t out < 10 GeV/c; 


8. 


Njet = 1; 




17. < 3%. 


* 


maximal P t of a hadron in the ECAL cell containing 


r a 7 dir -candidate; 




< A scalar sum of P t in the ring: P t sum (R = 0.4) - 


P t sum (i? = 0.2). 



Line number 1 of Table |7] makes primary preselection. It includes and specifies our first 
general cut (16) of Section 3.2 as well as the cut connected with ECAL geometry and the cut (19) 
that excludes r ) d%r -candidates accompanied by hadrons. 

Line number 2 of Table [7J fixes the values of Ptcur an d e cuT mat ' according to (17) and 
(18), define the isolation parameters of 7. 

The third cut selects the events containing 7 d * r -candidates with P t higher than ]3™ n (= 
CKIN(3)) threshold 46 . We impose the third cut to select the samples of events with Pt 1 > 40, 70 
and 100 GeV/c as ISR may smear the sharp kinematical cutoff defined by CKIN(3) ||9|- This 
cut reflects an experimental viewpoint when one is interested in how many events with -y dir - 
candidates are contained in some definite interval of P t 7 . 

The forth cut restricts a value of P^ g = P^ 0A ~ p ti?=o.2> where p tR l is a sum of P t of 
all ECAL cells contained in the cone of the radius R around the cell fired by 7* r -candidate ll34l . 

El. 

The fifth cut makes tighter the isolation criteria within R = 0.7 than those imposed onto 
7 djr -candidate in the second line of Table [7] 

The cuts considered up to now, apart from general preselection cut P t ^ et > 30 GeV/ c used 
in the first line of Table|7J were connected with photon selection ("photonic" cuts). Before we go 
further, some words of caution must be said here. Firstly, we want to emphasize that the starting 
numbers of the signal (S) and background (B) events (first line of Table EJ) may be specific only 
for PYTHIA generator and for the way of preparing primary samples of the signal and background 
events described above. So, we want to underline here that the starting values of S and B in the 
first columns of Table |8] are model dependent. 

But nevertheless, for our aim of investigation of efficiency of new cuts 11-17 (see fllO]- 
ifTTl ) the important thing here is that we can use these starting model numbers of S- and 5-events 
for studying a further relative influence of these cuts on S/B ratio. 

The cuts 6-9 are connected with the selection of events having only one jet and the defini- 
tion of jet-photon spatial orientation in 0-plane. The 9-th cut selects the events with jet and photon 

46 see PYTHIA manual 
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Fig. 13: Distribution of events over p t mlss in events with energetic e ±l s appearing as direct photon candidates for 
the cases Pt e > 70 GeV/c and Pt e > 100 GeV/c (here are used events satisfying cuts 1-5 of TableQ. 



transverse momenta being "back-to-back" to each other in 0— plane within the angle interval of 
the A0 = 17° size 47 . 

In line 10 we used the cut on p t miss to reduce a background contribution from the elec- 
troweak subprocesses qg — > q' + W ± and q q' — > g + with the subsequent decay W ± — > e ± u 
that leads to a substantial p f miss value. It is clear from the distributions over p t miss for two P t e 
intervals presented in Fig. [13] (compare with Fig. 13). One can see from the last column of Table 
|8]"e ± " that the cut on p t mtss reduces strongly (in about 4 times) the number of events containing 
e ± as direct photon candidate. 



Table 8: Values of significance and efficiencies for p J nm = 70 GeV/c. 



Cut 


S 


B 


Effs(%) 


Eff B (%) 


S/B 




1 


39340 


1247005 


100.00± 0.00 


100.000± 0.000 


0.03 


17562 


2 


36611 


51473 


93.06± 0.68 


4.128± 0.019 


0.71 


4402 


3 


29903 


18170 


76.01± 0.58 


1.457± 0.011 


1.65 


2038 


4 


26426 


11458 


67.17± 0.53 


0.919± 0.009 


2.31 


1736 


5 


23830 


7504 


60.57± 0.50 


0.602± 0.007 


3.18 


1568 


6 


23788 


7406 


60.47± 0.50 


0.594± 0.007 


3.21 


1554 


7 


23334 


6780 


59.3 1± 0.49 


0.544± 0.007 


3.44 


1460 


8 


19386 


4136 


49.28± 0.43 


0.332± 0.005 


4.69 


1142 


9 


18290 


3506 


46.49± 0.42 


0.281 ±0.005 


5.22 


796 


10 


18022 


3418 


45.81± 0.41 


0.274± 0.005 


5.27 


210 


11 


15812 


2600 


40.19± 0.38 


0.208± 0.004 


6.08 


176 


12 


13702 


1998 


34.83± 0.35 


0.160± 0.004 


6.86 


130 


13 


10724 


1328 


27.26± 0.30 


0.106± 0.003 


8.08 


88 


14 


10636 


1302 


27.04± 0.30 


0.104± 0.003 


8.17 


86 


15 


10240 


1230 


26.03± 0.29 


0.099± 0.003 


8.33 


84 


16 


8608 


984 


21.88± 0.26 


0.079± 0.003 


8.75 


64 


17 


6266 


622 


15.93± 0.22 


0.050± 0.002 


10.07 


52 



(*) The background (B) does not include the contribution from the "e events" (i.e. in which e ± fake 7-candidate) 
that is shown separately in the right-hand column "e*". 

Moving further we see from Table [U that the cuts 11-16 of Table [7] reduce the values of 
P t dust and P t out down to the values less than 10 GeV/c. The 17-th cut of Table □ imposes 

47 i.e. within the size of three calorimeter cells 



the jet isolation requirement. It leaves only the events with jets having the sum of P t in a ring 
surrounding a jet to be less than 3% of P t Jet . From comparison of the numbers in 10-th and 
17-th lines we make the important conclusion that all these new cuts (11-17), despite of model 
dependent nature of starting S/B value in line 10, may, in principle, lead to the following about 
two-fold improvement of S/B ratio. This improvement is reached by reducing the P t activity out 
of "7 + 1 jet" system. 

It is also rather interesting to mention that the total effect of "hadronic cuts" 6-17 for the 
case ofp™ m = 70 GeV/c consists of about twelve-fold decrease of background contribution at 
the cost of less than four-fold loss of signal events (what results in about 3.2 times growth of 
S/B ratio). So, in this sense, we may conclude that from the viewpoint of S/B ratio a study of 
"7 + jet" events may be more preferable as compared with a case of inclusive photon production. 

Table |9] includes the numbers of signal and background events left in three generated event 
samples after application of cuts 1-16 and 1-17. They are given for all three intervals of Pp '. 
Tables |9] and [8] are complementary to each other. The summary of Table [8] is presented in the 
middle section (pj nm = 70 GeV/c) of Table |9] where the line "Preselected" corresponds to the 
cut 1 of Tableland, respectively, to the line number 1 of Table[8]presented above. The line "After 
cuts" corresponds to the line 16 of Table [Hand line "+jet isolation" corresponds to the line 17 of 
Table El 



Table 9: Number of signal and background events remained after cuts. 



~ mm 
P± 

(GeV/c) 


Cuts 


7 

direct 


7 
brem 


photons from the mesons 




vr° 


V 


u 


A s 


40 


Preselected 


18056 


14466 


152927 


56379 


17292 


14318 


2890 


After cuts 


6238 


686 


824 


396 


112 


104 


24 


+ jet isol. 


3094 


264 


338 


150 


40 


44 


14 


70 


Preselected 


39340 


63982 


761926 


269666 


87932 


63499 


17562 


After cuts 


8608 


424 


320 


146 


58 


36 


64 


+ jet isol. 


6266 


262 


206 


90 


40 


24 


52 


100 


Preselected 


56764 


111512 


970710 


346349 


117816 


91416 


38872 


After cuts 


11452 


280 


124 


92 


24 


24 


136 


+ jet isol. 


9672 


204 


92 


64 


24 


20 


120 



Table 10: Efficiency, S/B ratio and significance values in the selected events without jet isolation cut. 





s 


B 


Eff S (%) 


Eff B (%) 


S/B 


s/Vb 


40 


6238 


2122 


34.55±0.51 


0.831±0.018 


2.9 


135.4 


70 


8608 


984 


21.88 ±0.26 


0.079 ± 0.003 


8.8 


274.4 


100 


11452 


544 


20.17 ±0.21 


0.033 ± 0.001 


21.1 


491.0 



Table 11: Efficiency, S/B ratio and significance values in the selected events with jet isolation cut. 



j5j™ n (GeV/c) 


s 


B 


Effs(%) 


Eff B (%) 


S/B 


s/Vb 


40 


3094 


836 


17.14±0.33 


0.327±0.011 


3.7 


107.0 


70 


6266 


622 


15.93 ± 0.22 


0.050 ± 0.002 


10.1 


251.2 


100 


9672 


404 


17.04 ±0.19 


0.025 ± 0.001 


23.9 


481.2 



Table[9]is done to show in more detail the origin of ^ dtr -candidates. The numbers in the "7— 
direct" column correspond to the respective numbers of signal events left in each of P// 1 intervals 



after application of the cuts defined in lines 1,16 and 17 of Table |7] (and column "5"' of Table 
lUl. Analogously the numbers in the "7 — brem" column of Table[9]correspond to the numbers of 
events with the photons radiated from quarks participating in the hard interactions. Columns 5-8 
of Table |9] illustrate the numbers of the "7 — mes" events with photons originating from ir°, 77, uj 
and Kg meson decays. In a case of Pf > 70 GeV/c the total numbers of background events, i.e. 
a sum over the numbers presented in columns 4-8 of Tabled are shown in the lines 1, 16 and 17 
of column "5" of Table El The other lines of Table for p^ in = 40 and 100 GeV/c have the 
meaning analogous to that described above for p j™ n = 70 GeV j c. 

The last column of Table[9] shows the number of preselected events with e^. 

The numbers in Tables flOl (without jet isolation cut) and [TT] (with jet isolation cut) accumu- 
late in a compact form the final information of Tables |7]-|9l Thus, for example, the columns S 
and B of the line that corresponds to p ™ m = 70 GeV/ c contain the total numbers of the selected 
signal and background events taken at the level of 16-th (for Table [KB and 17-th (for Table fill 
cuts from Table 

It is seen from Table flOl that in the case of Selection 1 the ratio S/ B grows from 2.9 to 21.1 
while P? increases from P$ > 40 GeV/c to P? > 100 GeV/c interval. 

The jet isolation requirement (cut 17 from Table|71) noticeably improves the situation at low 
Pt' (see Table fTTTl. After application of this criterion the value of S/B increases from 2.9 to 3.7 
at Pp > 40 GeV jc and only from 21.1 to 23.9 at Pt' > 100 GeV jc. Remember on this occasion 
the conclusion that the sample of events selected with our criteria has a tendency to contain more 
events with an isolated jet as P t 7 increases 48 . 

Let us underline here that, in contrast to other types of background, "7 — brem" background 
has an irreducible nature. Thus, the number of "7 — brem" events should be carefully estimated 
for each Pt' interval using the particle level of simulation in the framework of event generator 
like PYTHIA. They are also have to be taken into account in experimental analysis of the prompt 
photon production data at high energies. 

Tables fT2l and IT3l shows the relative contributions of fundamental QCD subprocesses (hav- 
ing the largest cross sections) qg — > qg, qq — > qq, gg — > qq and gg — > gg 49 that define the main 
production of "7— brem" and '7— mes" background in event samples selected with criteria 1-13 
of Table |7]in three Pt 1 intervals. 

Accepting the results of simulation with PYTHIA, we found from the event listing analysis 
that in the main part of selected "7 — brem" events these photons are produced in the final state 
of the fundamental 2 — > 2 subprocess 5() . Namely, they are mostly radiated from the outgoing 
quarks in the case of the first three sets of subprocesses {qg — > qg, qq —> qq and gg — > qq). They 
may also appear as a result of string breaking in a final state of gg — > gg scattering. But this 
subprocess, naturally, gives a small contribution into "7 + jet" events production. 

It may be noted also from the first two columns of Tables fT2l and fT~3l that the most of "7 — 
brem" and "7— mes" background events (93% at least) originate from qg — > qg and g^- — > q^qj, 
QiQi — * qj<lj subprocesses. Tables fT2l and fT~3l show also a tendency of increasing the contribution 
from the subprocess qiqj q^qj and q^ — > q^ (given in the second columns of tables) with 
growing P t 7 . 

The values of signal-to-background ratios in Tables ITOl ITTI are obtained without any detec- 
tor effects. But these numbers can be noticeably increased if we take into account information 

48 see Sections 5-7 and Appendix 2 
49 ISUB=1 1, 12, 28, 53 and 68 (see 0) 
50 i.e. from lines 7, 8 in Fig. 3 



Table 12: Relative contribution (in per cents) of different QCD subprocesses into the "7 — brem" events production. 



Pi' 

{GeV/c) 


fundamental QCD subprocess 


Q9 Q9 


qq -> gg 


55 -> <?g 


55 -> 55 


40-70 


62.1±6.6 


31.8±4.0 


3.3±1.0 


2.8±0.9 


70-100 


52.3±7.7 


42.4±6.4 


3.8±1.4 


1.5±0.9 


> 100 


41.8±6.0 


56.9±7.2 


1.3±0.7 





Table 13: Relative contribution (in per cents) of different QCD subprocesses into the "7— mes" events production. 



(GeV/c) 


fundamental QCD subprocess 


95 -> 95 


gg -> gg 


55 -> 99 


55 -> 55 


40-70 


59.3±5.2 


34.8±3.5 


2.9±0.7 


2.4±0.7 


70-100 


48.6±8.0 


47.3±7.8 


0.5±0.5 


0.5±0.5 


> 100 


41.8±6.4 


53.9±7.6 


1.8±0.9 


0.7±0.5 



from the preshower detector 51 . First of all, photons in the signal "7 + jet" events have the distri- 
bution over number of preshower 3-dimensional 52 clusters N^ st different from one for photon 
candidates in the QCD background events. Selection efficiencies for the signal and background 
events after application of the cut N^ st < 1 is shown in Fig. 1141 for \r]^\ < 0.9. Relatively big 
numbers of N^ st in the QCD background may be explained by the facts that besides 7r°'s we 
have a contribution from events with multiphoton decays of 77, K® and u mesons and that despite 
the strong isolation criteria photon candidates from the background events still have a hadron 
accompaniment. 

Additional rejection can be obtained after analysis of energy distributions among the strips 
of each of three single layer clusters (SLC). They are again different for the signal and back- 
ground events. As parameters for the discrimination one can take the energy weighted widths of 
three SLC's and ratios of energy deposited in the hottest strip to the total energy of SLC cluster 
E™* I Eslc- The selection efficiencies for single 7's and 7r°'s (as a most difficult case from the 
point of view of discrimination) are presented in Fig. [131 53 . 

Thus, the total effect of data analysis in the preshower detector can lead to additional in- 
crease in the S/ B of order of 3 — 4 54 . 

From Tables l9l - ITT1 we have seen that the cuts listed in Table |7] (having rather moderate 
values of Ptcur an d Ptcur) allow to suppress the major part of the background events. The 
influence of these two cuts on: 

(a) the number of selected events (for L int = 300 pb^ 1 ); 

(b) the signal-to-background ratio S/B; 

(c) the mean value of F = [Pi —P t Jet ) / 'Pt and its standard deviation value a(F) 
is presented in Tables 1-12 of Appendix 3 for their variation in a wide range. 

Let us emphasize that the tables of Appendix 3 include, in contrast to Appendix 2, the re- 

51 Central (CPS) and forward (FPS) preshower detectors are placed at [77] < 1.1 and 1.2 < |?/| < 2.5, respectively, 
and have a similar 3-layered architecture with set of triangular scintillator strips in every layer. 

52 because they are built from 3 layers rotated in the space by some angles with respect to each other 

53 A consideration of the full QCD background left after our selection cuts (see cuts 1 — 16 of Table 0plus re- 
quirement N^ugt < 1 above) is very difficult because of a pure statistics. Obviously should decrease "7 — mes" 
background selection due to the "77, K®,u> events" contribution and probably due to an admixture of hadron accom- 
paniment around 7-candidates in those events. 

54 these factors are caused mainly by the single photon selection efficiency and P t 7 interval. 
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Fig. 14: Selection efficiencies for photons from "7 + jet" process and photon candidates from QCD background 
obtained after cut on the number of 2-D clusters in the central preshower: N^ st < 1. 
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Fig. 15: Selection efficiency of single photon via selection efficiency of 7r° obtained by using two set of quantities, 
measured in the preshower detector: three energy weighted widths of the single layer clusters (full line) and the same 
plus three ratios of energy deposited in the hottest strip to the total energy of SLC clusters I Es lc (dashed 
line). Four P t values, 20, 40, 60, 80, are considered on the plots above. 



suits obtained after analyzing three generated samples (described in the beginning of this section) 
of signal and background events. These events were selected with the cuts of TableEJ 

Namely, the cuts 1-10 of Table|7]were applied for preselection of "7 + 1 jet" events. The 
jets in these events as well as clusters were found by use of only one jetfinder LUCELL (for the 
whole 77 region \r]^ et \ < 4.2). 

Tables 1-4 of Appendix 3 correspond to the simulation with pj™ n = 40 GeV/c. Anal- 
ogously, the values of p ™ m = 70 GeV / c and p J nm = 100 GeV / c were used for Tables 5-8 
and Tables 9-12 respectively. The rows and columns of Tables 1-12 illustrate, respectively, the 
influence of Ptcur an ^ Ptcur on me quantities mentioned above in the points (a), (b), (c). 

First of all, we see from Tables 2, 6 and 10 of Appendix 3 that a noticeable reduction of 
the background take place while moving along the table diagonal from the right-hand bottom 
corner to the left-hand upper one, i.e. with reinforcing Ptcur an d Ptcur- So, we see that for 
pmm _ 4Q Q e y/ C the value of S/B ratio changes in the table cells along the diagonal from 
S/B = 2.2 (in the case of no limits on these two variables), to S/B = 2.9 for the cell with 
p t <gff* = 10 GeV/c and Ptcur = 10 GeV/c. Analogously, for p^ in = 100 GeV/c the value of 
S/ B changes in the same table cells from 10.0 to 29.5 (see Table 10 of Appendix 3) 55 . 

The second observation from Appendix 3. The restriction of Ptcur an d Ptcur improves 
the calibration accuracy. Table 3 shows that in the interval P^ > 40 GeV/ c the mean value of the 
fraction F(= (Pt'-P t Jet ) / Pt') decreases from 0.049 (the bottom right-hand corner) to 0.024 for 
the table cell with P t cur = 10 GeV/c and P t °Sur = 10 GeV/c. At the same time, the both cuts 
lead to a noticeable decrease of the Gaussian width cr(F) (see Table 4 and also Tables 8, 12). For 
instance, for p™ m = 40 GeV /c a(F) drops by about a factor of two: from 0.159 to 0.080. It 
should be also noted that Tables 4, 8 and 12 demonstrate that for any fixed value of Ptcur further 
improvement in a(F) can be achieved by limiting P t out (e.g. in line with P t $tf* = 10 GeV/c 
a(F) drops by a factor of 2 with variation of P t out from 1000 to 5 GeV/c). 

The explanation is simple. The balance equation (l25t contains 2 terms on the right-hand 
side (1 — cosAcp) and P t (0+r) > 4.2)/P t 7 . The first one is negligibly small in a case of Selection 
1 and tends to decrease with growing P t 7 (see tables in Appendix 2). So, we see that in this case 
the main source of the disbalance in equation (l25t is the term P t (0+7] > 4.2)/ P^ . This term can 
be diminished by decreasing P t activity beyond the jet, i.e. by decreasing P t 



out 



The behavior of the number of selected events (for L int = 300 pb x ), the mean values 
of F = (PS - P t Jet ) I Pt' and its standard deviation a(F) as a function of P t ° c u UT (with fixed 
Ptcur = 10 GeV /c) are also displayed in Fig.[l6]for events with non-isolated (left-hand column) 
and isolated jets (right-hand column, see also Tables 13-24 of Appendix 3). 

Thus, we can conclude that application of two criteria introduced in Section 3.2, i.e. Ptcur 
and Pt°cur> results in two important consequences: significant background reduction and essen- 
tial improvement of the calibration accuracy. 

The numbers of events (for L int = 300 pb" 1 ) for different Ptcur an ^ Ptcur are gi ven m 
the cells of Tables 1, 5 and 9 of Appendix 3. One can see that even with such strict Ptcur an ^ 
Ptcur values as, for example, 10 GeV/c for both we would have a sufficient number of events 
(about 100 000, 7 000 and 1 300 for P? > 40 GeV/c, P{< > 70 GeV/c and > 100 GeV/c, 
respectively) with low background contamination (S/B = 2.9, 8.8 and 21.1) and a good accuracy 
of the PS - P t Jet balance: F = 2.4%, 1.5% and 1.2%, respectively, for the case of Selection 1. 

In addition, we also present Tables 13-24 of Appendix 3. They contain the information 

55 even better results produces a combined application of stronger cuts Ptcur = 5 GeV/ c and Ptcur = 5 GeV/c 
(see Appendix 3) 
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Fig. 16: Number of events (for Lj„i = 300 pb^ 1 ), mean value of (Pf 1 -P t Jet )/Pi l (= F) and its standard deviation 
u{F) distributions over p t out for the cases of non-isolated (left-hand column) and isolated (right-hand column) jet 
and for three intervals: P t *> > 40, 70 and 100 GeV/c. PtcuT = 10 GeV/c. 



analogous to that in Tables 1-12 but for the case of isolated jets with e je * < 3%. From these 
tables we see that with the same cuts PtcvT = ^Wut = 10 GeV/c one can expect about 
47 000, 5 000 and 1000 events for > 40 GeV/c, > 70 GeV/c and > 100 GeV/c, 
respectively, with a much more better fractional PP - P t Jet balance: F = 0.5%, 0.7% and 0.1%. 

Let us mention that all these PYTHIA results give us an indication of a tendency and may 
serve as a guideline for further full GEANT simulation that would allow to come to a final con- 
clusion. 

To conclude this section we would like to stress, firstly, that, as is seen from Tables the 
"7 — brem" background defines a dominant part of the total background. One can see from Table 
|9]fhat 7T° contribution being about the same as "7 — brem" at p J nm > 40 GeV/c becomes three 
times less than "7 — brem" contribution at p™ m > 100 GeV/c. We would like to emphasize here 
that this is a strong prediction of the PYTHIA generator that has to be compared with predictions 
of another generator like HERWIG, for example. 

Secondly, we would like to mention also that, as it is seen from Tables [8] and the photon 
isolation and selection cuts 1-6, usually used in the study of inclusive photon production (see, 
for instance, 01, 01, 03), increase the S/B ratio up to 3.20 only (for > 70 GeV/c). The 



other cuts 6-17, that select events with a clear "7 + jet" topology and limited P t activity beyond 
"7 + jet" system, lead to quite a significant improvement of S/B ratio by a factor of three (to 
S/B = 10.07). 

The numbers in the tables of Appendix 3 were obtained with inclusion of the contribu- 
tion from the background events. The tables show that their account does not spoil the PP — 
P t i et balance in the event samples preselected with the cuts 1-10 of Table|7| The estimation of the 
number of these background events would be important for the gluon distribution determination 
(see Section 9). 

9. "7 + jet" EVENT RATE ESTIMATION FOR GLUON DISTRIBUTION DETERMI- 
NATION AT THE TEVATRON RUN II. 

The number of "7 + jet" events suitable for measurement of gluon distribution in different x and 
Q 2 intervals at Run II is estimated. It is shown that with Li n t = 3 fb^ 1 it would be possible to collect 
about one million of these events. This number would allow to cover a new kinematical area not studied in 
any previous experiment (10~ 3 < x < 1.0 with 1.6 • 10 3 < Q 2 < 2 • 10 4 (GeV/c) 2 ). This area in the region 
of small x > 10~ 3 has Q 2 by about one order of magnitude higher than reached at HERA now. 

As many of theoretical predictions for production of new particles (Higgs, SUSY) at the 
Tevatron are based on model estimations of the gluon density behavior at low x and high Q 2 , the 
measurement of the proton gluon density for this kinematic region directly in Tevatron experi- 
ments would be obviously useful. One of the promising channels for this measurement, as was 
shown in QUI , is a high P t direct photon production pp{p) — > r ) dlT + X. The region of high 
P u reached by UA1 ED, UA2 01, CDF E3 and DO El extends up to P t « 80 GeV/c and 
recently up to P t = 105 GeV/c 051 . These data together with the later ones (see references in 
ll37l - ll45l and recent E706 ll46l and UA6 ll47ll results) give an opportunity for tuning the form of 
gluon distribution (see EH)- The rates and estimated cross sections of inclusive direct 

photon production at the LHC were given in OOl (see also [[49]). 

Here for the same aim we shall consider the process pp — > r y d%r + 1 Jet + X defined in the 
leading order by two QCD subprocesses (la) and (lb) (for experimental results see l!50l . iBTlO . 

Apart from the advantages, discussed in Section 8 in connection with the background sup- 
pression (see also tl52l-[l58l ). the "7^ + 1 Jet" final state may be easier for physical analysis 
than inclusive photon production process "7 dir + X" if we shall look at this problem from the 
viewpoint of extraction of information on the gluon distribution in a proton. Indeed, in the case 
of inclusive direct photon production the cross section is given as an integral over the products 
of a fundamental 2^2 parton subprocess cross sections and the corresponding parton distri- 
bution functions f a (x a , Q 2 ) (a = quark or gluon), while in the case of pp — > i y dir + 1 Jet + X 
for P t Jet > 30 GeV/c (i.e. in the region where "k t smearing effects" 56 are not important, see 
ll43l ) the cross section is expressed directly in terms of these distributions (see, for example, |41 ]): 

, -p- , ~ 2 = x afa(x a ,Q 2 ) x b f b (x b ,Q 2 )^(ab^ cd), (30) 
d^d^dPt ab dt 

where 

x a ,b = Pt/Vs- (exp(±r] 1 ) + exp(±7] 2 )). (31) 

The designation used above are as the following: rj\ = rp , r]2 = r] Jet ; P t = Pt\ a,b = 
q, q, g; c,d = q, q, g, 7. Formula (l30l and the knowledge of q, q distributions allow the gluon dis- 

56 This terminology is different from ours, used in Sections 2 and 9, as we denote by "fc t " only the value of parton 
intrinsic transverse momentum. 



tribution f g (x, Q 2 ) to be determined after account of selection efficiencies for jets and ^ d ' lT — candidates 
as well as after subtraction of the background contribution, left after the used selection cuts 1-13 
of Table|7](as it was discussed in Section 8 keeping in hand this physical application). 

In the previous sections a lot of details connected with the structure and topology of these 
events and the features of objects appearing in them were discussed. Now with this information 
in mind we are in position to discuss an application of the "7 + jet" event samples, selected 
with the previously proposed cuts, for estimating the rates of the gluon-based subprocess (la) in 
different x and Q 2 intervals. 

Table PHI shows percentage of "Compton-like" subprocess (la) (amounting to 100% to- 
gether with (lb)) in the samples of events selected with cuts (16)-(22) of Section 3.2 for Ptcur = 
10 GeV/c for different P t 7 and r] jet intervals: Central (CC) (\r] jet \ < 0.7) 57 , Intercryostat (IC) 
0.7 < \7]^ et \ < 1.8 and End (EC) 1.8 < |?7 je '| < 2.5 parts of calorimeter. We see that the contribution 
of Compton-like subprocess grows by about 5 — 6% with \rj^ et \ enlarging and drops with growing 
P t i et (tt P t 7 in the sample of the events collected with the cuts 1 — 13 of Table ITJ . 



Table 14: The percentage of Compton-like process q g — ► 7 + q. 



Calorimeter 
part 


P t Jet interval (GeV/c) 


40-50 


50-70 


70-90 


90-140 


CC 


84 


80 


74 


68 


IC 


85 


82 


76 


70 


EC 


89 


85 


82 


73 



In Table IT31 we present distribution of the number of events that are caused by the q g — > 
7 + q subprocess, in various intervals of the Q 2 (= (Pi 1 ) 2 ) 58 and x (defined according to (OTT) '). 
These events have passed the following cuts (P t out was not limited): 

P 7 > 40 GeV/c, \if \ < 2.5, P t Jet > 30 GeV/c, \ V Jet \ < 4.2, P t hadr > 7 GeV/c, 

P^ut = 4 GeV/c, el UT = 7%, A0 < 17°, = 10 GeV/c. (32) 



Table 15: Number of g q — > 7™" + q events at different Q 2 and x intervals for Li nt =3/6 



Q 2 

(GeV/c) 2 


x values of a parton 


All x 


.001 - .005 


.005 - .01 


.01 - .05 


.05 -.1 


.1 - .5 


.5 - 1. 


.001 - 1. 


1600-2500 


8582 


56288 


245157 


115870 


203018 


3647 


632563 


2500-4900 


371 


13514 


119305 


64412 


119889 


3196 


320688 


4900-8100 





204 


17865 


13514 


26364 


1059 


59007 


8100-19600 








3838 


5623 


11539 


548 


21549 




1 033 807 



The analogous information for events with the charmed quarks in the initial state g c —> 
^dir _|_ c - IS p resen t e d i n Table fT6l The simulation of the process g b — > ^ dlT + b has shown that 
the rates for the 6-quark are 8 - 10 times smaller than for the c-quark. These event rates are also 
given in Appendix 1 for different P^ intervals in the lines denoted by "Neventsfc/b)" 59 . 

57 see also tables of Appendix 1 containing lines "29sub/all" 
58 see 

59 Analogous estimation for LHC energy was done in ITHl and (59). 



Table 16: Number of g c — > ^ dlr + c events at different Q 2 and x intervals for Li nt = 3/6 1 . 



Q 2 
{GeV/cf 


x values of a parton 


Alls 


.001 - .005 


.005 - .01 


.01 - .05 


.05 -.1 


.1 - .5 


.5-1. 


.001 - 1. 


1600-2500 


264 


2318 


21236 


11758 


14172 


58 


49805 


2500-4900 


13 


332 


9522 


6193 


7785 


40 


23885 


4900-8100 





4 


914 


1055 


1648 


16 


3637 


8100-19600 








142 


329 


612 


8 


1092 




78 419 



Fig. 17 shows in the widely used 
(x, Q 2 ) kinematic plot (see ll60ll and also 
in 114310 what area can be covered by 
studying the process q g — > 7 + q at 
Tevatron. The distribution of number 
of events in this area is given by Ta- 
ble [15] From this figure and Table fT5l it 
becomes clear that with integrated lumi- 
nosity L int = 3 fb" 1 it would be possi- 
ble to study the gluon distribution with 
a good statistics of "7 + jet" events in 
the region of 10~ 3 < x < 1.0 with Q 2 
by about one order of magnitude higher 
than reached at HERA now. It is worth 
emphasizing that extension of the exper- 
imentally reachable region at the Teva- 
tron to the region of lower Q 2 overlap- 
ping with the area covered by HERA 
would also be of great interest. 

Fig. 17: The (x, Q 2 ) kinematic region for studying pp — > 7 + Jet process at Tevatron Run II. 



10. SUMMARY. 

We have done an attempt here to consider, following lfTUll-[|T8l. the physics of high P t direct 
photon and jet associative production in proton-antiproton collisions basing on the predictions of 
PYTHIA generator and the models implemented there. This work may be useful for two practical 
goals: for absolute jet energy scale determination and for gluon distribution measurement at 
Tevatron energy. 

The detailed information provided in the PYTHIA event listings allows to track the origin of 
different particles (like photons) and of objects (like clusters and jets) that appear in the final state. 
So, the aims of this work was to explore at the particle level as much as possible this information 
for finding out what effect may be produced by new variables, proposed in H1OI - 017I for selection 
of "7 + jet" events, and the cuts on them for solution of the mentioned above practical tasks. 




For the first problem of the jet energy determination an important task is to select the events 
that may be caused (with a high probability) by the qq — > g + 7 and qg — > q + 7 fundamental 
parton subprocesses of direct photon production. To take into account a possible effect of initial 
state radiation (its spectra are presented in different P 4 7 intervals in Section 5) we used here 
the P t -balance equation (see (fT5t ) written for an event as a whole. It allows to express — 
P t i et fractional disbalance (see (12*51) ) through new variables lfTDl -lfT7l that describe the P t activity 
out of "7 + jet" system. They are P t out and p t clus \ i.e. P t of mini-jets or clusters that are 
additional to the main jet in event. The latter is the most "visible" part of P° ut . 

The sources of the P t 7 — P/ e * disbalance are investigated. It is shown that the limitation of 
P t of clusters, i.e. p t dust , can help to decrease this disbalance. Analogously, the limitation of P t 
activity of all detectable particles (\r]i\ < 4.2) beyond the "7 + jet" system, i.e. P t out , also leads 
to a noticeable P t 7 — P t jet disbalance reduction (see Sections 7,8). 

It is demonstrated that in the events selected by means of simultaneous restriction from 
above of the p t dust and P t out activity the values of Pf and P t jet are well balanced with each 
other. The samples of these "7 + jet" events gained in this way are of a large enough volume 
for the jet energy scale determination in the interval 40 < P t 7 < 140 GeV/c (see Tables 1-12 of 
Appendix 3). 

It is worth mentioning that the most effect for improvement of P t 7 and P^ et balance can 
be reached by applying additionally the jet isolation criterion defined in [fTOj — [|T7l . As it can be 
seen from Tables 13-18 of Appendix 2 and Tables 13-24 of Appendix 3, the application of this 
criterion allows to select the events having the P{' — P t jet disbalance at the particle level less 
than 1%. Definitely, the detector effects may worsen the balance determination due to a limited 
accuracy of the experimental measurement 60 . 

We present also PYTHIA predictions for the dependence of the distributions of the number 
of selected "7 + jet" events on P t 7 and rf^ (see Section 5 and also tables of Appendix 2 with 
account of P t dust variation). 

The corrections to a jet P t the measurable values of P^ et that have take into account the 
contribution from neutrinos belonging to a jet are presented for different P t Je *(^ P t 7 for the 
selected events) intervals in the tables of Appendix 1. It is shown in Section 4 that a cut on 
pmiss < Q e yj c allows to reduce the neutrino contribution to the value of (Ptf^) a u events = 
0.1 GeV/c. 

At the same time, as it is shown in ll2"6*1. and discussed in Section 8 (see also JT71 ). this 
cut noticeably decreases the number of the background e^ -events in which e ± (produced in the 
W ± — > e ± v weak decay) may be registered as direct photon. 

The possibility of the background events (caused by QCD subprocesses of qg, gg, qq scat- 
tering) suppression was studied in Section 8. Basing on the introduced selection criteria that 
include 17 cuts (see Table [7] of Section 8), the background suppression relative factors and the 
values of signal event selection efficiencies are estimated (see Tables [8lfTTb. 

It is shown that after applying the first 6 "photonic" cuts (that may be used, for example, 
for selecting events with inclusive photon production and lead to S/B ratio equal to 3.2 for 
Pt 1 > 70 GeV/ c, see Table [U) the use of the next 1 1 "hadronic" cuts of Table [7] may lead to 
further essential improvement of Sj B ratio (by factor of 3.2 for the same P t 7 > 70 GeV j c where 
S/B becomes 10.1. 

60 We are planning to present the results of full GEANT simulation with the following digitization and reconstruc- 
tion of signals by using the corresponding DO packages in the forthcoming papers. 



It is important to underline that this improvement is achieved by applying "hadronic" cuts 
that select the events having clear "7 + jet" topology at the particle level and also having 
rather "clean" area (in a sense of limited P t activity) beyond a "7 + jet" system. In this sense 
and taking into account the fact that these "hadronic" cuts lead to an essential improvement of 
PP — P^ et balance, one may say that the cuts on p t dust and P° ut , considered here, do act quite 
effectively to select the events caused by leading order diagrams (see Fig. 1) and do suppress the 
contribution of NLO diagrams, presented in Figs. 2, 4. 

The consideration of the cuts, connected with detector effects (e.g., based the preshower 
usage) may lead to further noticeable improvement of S/B ratio. 

Another interesting predictions of PYTHIA is about the dominant contribution of "7-brem" 
events into the total background at Tevatron energy, as in was already mentioned in Section 8 (see 
also ifPTI and ll26l ). As the "7-brem" background has irreducible nature its careful estimation is 
an important task and requires the analogous estimation with another generator. 

To finish the discussion of the jet calibration study let us mention that the main results on 
this subject are summed up in Tables 1-12 (Selection 1) and 13-24 (Selection 2 with jet isolation 
criterion) of Appendix 3 and Fig. fT*6l 

It should be emphasized that numbers presented in all mentioned tables and figures were 
found within the PYTHIA particle level of simulation. They may depend on the used generator 
and on the particular choice of a long set of its parameters 61 as well as they may change after 
account of the results of the full GEANT-based simulation. 

It is also shown that the samples of the "7 + jet" events, gained with the cuts used for 
the jet energy calibration, can provide an information suitable also for determining the gluon 
distribution inside a proton in the kinematic region (see Fig. 17) that includes x values as small as 
accessible at HERA iRHl. [R>2"1. but at much higher Q 2 values (by about one order of magnitude): 
1(T 3 < x < 1.0 with 1.6 • 10 3 < Q 2 < 2 • 10 4 (GeV/c) 2 . The number of events, based on the 
gluonic process (la), that may be collected with L int = 3 fb^ 1 in different x- and Q 2 - intervals 
of this new kinematic region for this goal are presented in Table [T31 (all quarks included) and in 
Table fT6l (only for charm quarks). 
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We have already mentioned that we are planning to perform analogous analysis by help of another generator like 
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Appendix 1 

Selection I. </) (7jet) = 180° ±17°. UAl algorithm. L int = 300 pb 



Table 1: 40 < PC < 50 GeV/c. 



r> clust 
r tcUT 


30 


20 


15 


10 


5 


p t jet 


42.646 


42.460 


42.410 


42.564 


42.912 


P t JGt _P t J et 


0.127 


0.127 


0.131 


0.133 


0.105 


i) Jet 


0.129 


0.128 


0.133 


0.135 


0.106 


TiVdJet 

-^event 


0.171 


0.170 


0.169 


0.166 


0.152 


j ) Jet 


0.153 


0.153 


0.158 


0.157 


0.113 


event 


0.148 


0.146 


0.146 


0.144 


0.126 


t-> miss 
"t 


2.088 


2.083 


2.096 


2.105 


2.101 


p miss 


2.366 


2.370 


2.383 


2.403 


2.310 


Nevent( c ) 


964 


926 


865 


723 


348 


Nevent^) 


100 


94 


90 


70 


34 


29sub/all 


0.85 


0.85 


0.84 


0.84 


0.83 


Entries 


10493 


10144 


9472 


7992 


4421 



Table 2: 70 < PC < 90 GeV/c. 



p clust 
r tcUT 


30 


20 


15 


10 


5 


p t jet 


72.873 


74.375 


75.239 


75.968 


76.353 


p Jet_p jet 


0.257 


0.259 


0.272 


0.250 


0.245 


\) Jet 


0.259 


0.262 


0.275 


0.253 


0.248 


r>v£Jet 
-^event 


0.182 


0.176 


0.177 


0.175 


0.173 


tj Jet 


0.184 


0.181 


0.186 


0.168 


0.174 


TjHEJet 
-^event 


0.172 


0.169 


0.171 


0.172 


0.165 


t-> rmss 
"t 


2.178 


2.182 


2.196 


2.168 


2.190 


p miss 


3.092 


3.123 


3.179 


3.118 


3.089 


Nevent( c ) 


129 


108 


91 


64 


30 


Nevent^ 


22 


18 


13 


9 


2 


29sub/all 


0.77 


0.76 


0.76 


0.75 


0.72 


Entries 


13641 


11613 


9892 


7495 


3845 



Table 3: 90 < P t 7 < 140 GeV/c. 



r> clust 
r tCUT 


30 


20 


15 


10 


5 


P/ et 


101.878 


103.159 


103.988 


104.565 


104.615 


p Jet_p jet 


0.331 


0.330 


0.319 


0.312 


0.317 


1 ) Jet 
Pt(u) 


0.334 


0.333 


0.321 


0.315 


0.320 


TjvGJet 
-^event 


0.190 


0.188 


0.187 


0.185 


0.179 


I ) Jet 

F Hn) 


0.272 


0.283 


0.272 


0.280 


0.309 


TjuEJet 
event 


0.181 


0.180 


0.175 


0.170 


0.163 


1 1 miss 
"t 


2.186 


2.197 


2.193 


2.195 


2.201 


p miss 
r t U £jet 


3.339 


3.339 


3.276 


3.238 


3.345 


Nevent( e ) 


51 


40 


32 


22 


9 


Nevent^) 


6 


5 


4 


2 


1 


29sub/all 


0.70 


0.70 


0.69 


0.68 


0.66 


Entries 


14058 


11806 


9997 


7439 


3673 



Appendix 2 P t lso1 < 4 GeV/c, e< < 7%, 

(j>(yjet) = 180° ± 17°. UA1 algorithm. L int = 300 pb' 1 . 



Table 4: Selection 1. 40 < P 7 < 50 GeV/c. 



r> clust 
r tcUT 


30 


20 


15 


10 


5 


Nevent 


12915 


12486 


11659 


9837 


5442 


P56 


10.1 


9.6 


8.9 


7.9 


6.2 


A<p 


5.9 


5.8 


5.6 


5.1 


3.9 


p out 


7.6 


7.3 


6.8 


6.1 


4.6 


p |T7| >4.2 


2.0 


2.0 


2.0 


2.0 


2.0 


(Pj-P t P a «)/PJ 


0.0120 


0.0155 


0.0147 


0.0116 


0.0071 


(P t J -Pt Part )/P t J 


-0.0291 


-0.0291 


-0.0296 


-0.0275 


-0.0213 


(Pt 1 -Pt J )/Pt 7 


0.0363 


0.0400 


0.0400 


0.0357 


0.0266 


P(0+77>5)/P t 7 


0.0279 


0.0319 


0.0325 


0.0293 


0.0226 


1 - cos(A0) 


0.0084 


0.0081 


0.0076 


0.0064 


0.0040 


<7(£>&l7,JJ) 


0.1531 


0.1414 


0.1298 


0.1142 


0.0904 


cr(£>6[7,part]) 


0.1957 


0.1831 


0.1667 


0.1424 


0.1105 


Entries 


10493 


10144 


9472 


7992 


4421 



Table 5: Selection 1. 70 < P 7 < 90 GeV/c. 



r> clust 
r tcUT 


30 


20 


15 


10 


5 


Nevent 


2414 


2055 


1751 


1327 


681 


P56 


14.7 


12.5 


11.0 


9.1 


6.8 


A<f> 


5.4 


4.7 


4.2 


3.4 


2.3 


p out 


12.5 


10.4 


8.9 


7.0 


4.9 


p t M>4.2 


1.9 


1.9 


2.0 


2.0 


2.0 


(P t T-P t P°")/P t T 


0.0328 


0.0184 


0.0118 


0.0067 


0.0038 


(Pt J -P t part )/P t J 


-0.0411 


-0.0310 


-0.0244 


-0.0192 


-0.0151 


(Pt 1 ~Pt J )/Pt J 


0.0642 


0.0440 


0.0325 


0.0233 


0.0171 


P t (0+?7>5)/P t 7 


0.0570 


0.0382 


0.0279 


0.0203 


0.0156 


1 - cos(A4>) 


0.0073 


0.0058 


0.0046 


0.0030 


0.0014 


a(Db[ 1 ,J\) 


0.1518 


0.1207 


0.1015 


0.0812 


0.0624 


a(Db[y,part) 


0.1789 


0.1467 


0.1268 


0.1048 


0.0843 


Entries 


13641 


11613 


9892 


7495 


3845 



Table 6: Selection 1. 90 < P t 7 < 140 GeV/c. 



r> clust 
^tCUT 


30 


20 


15 


10 


5 


Nevent 


1242 


1043 


885 


669 


333 


P56 


15.0 


12.7 


11.2 


9.4 


7.0 


Acf) 


4.5 


3.7 


3.2 


2.5 


1.8 


D out 


13.2 


10.6 


9.0 


7.1 


5.0 


p t M>4.2 


1.9 


1.9 


1.9 


1.9 


1.9 


(Pt~<-Pt Part )lPt~< 


0.0102 


0.0045 


0.0014 


0.0007 


0.0003 


{ p.J_p r rt )/p .J 


-0.0382 


-0.0276 


-0.0221 


-0.0160 


-0.0121 


(P 7 -P J )/P 7 


0.0417 


0.0286 


0.0213 


0.0153 


0.0112 


P t {0+ V >5)/Pt J 


0.0363 


0.0248 


0.0185 


0.0136 


0.0103 


1 - cos(A(f>) 


0.0054 


0.0038 


0.0028 


0.0018 


0.0009 


a(Db[ 7 ,J\) 


0.1154 


0.0896 


0.0753 


0.0605 


0.0479 


a(Db[y,part]) 


0.1359 


0.1111 


0.0981 


0.0861 


0.0677 


Entries 


26759 


22471 


19068 


14411 


7163 



'Number of events (Nevent) is given in this and in the following tables for integrated luminosity L; nt = 300 pb 1 . 



Table 7: Selection 2. 40 < P t 7 < 50 GeV/c (e^ 4 < 3%. 



^tCUT 


30 


20 


15 


10 


5 


Nevcnt 


5189 


5043 


4804 


4222 


2689 


P56 


9.4 


8.9 


8.4 


7.4 


5.9 


At/> 


5.5 


5.4 


5.3 


4.8 


3.8 


/ > out 


7.0 


6.7 


6.3 


5.6 


4.3 


p t \n\>4.2 


2.0 


2.0 


2.0 


2.0 


2.0 


{p t -<-p t p an )/p t -< 


-0.0237 


-0.0179 


-0.0143 


-0.0126 


-0.0085 


(Pt j -Pt part )/Pt j 


-0.0078 


-0.0094 


-0.0105 


-0.0135 


-0.0125 


{Pt 1 -Pt^/Pt 1 


-0.0163 


-0.0088 


-0.0043 


0.0001 


0.0032 


p t (o+n>5)/p t - y 


-0.0238 


-0.0161 


-0.0111 


-0.0058 


-0.0005 


1 - cos(A(f>) 


0.0076 


0.0074 


0.0069 


0.0059 


0.0038 


<T(Db[ 7 ,J\) 


0.1531 


0.1373 


0.1253 


0.1082 


0.0878 


a{Db[^i,part\) 


0.1814 


0.1661 


0.1515 


0.1251 


0.1028 


Entries 


4216 


4097 


3903 


3430 


2185 



Table 8: Selection 2. 70 < P t 7 < 90 GeV/c (e jet < 3%. 



r> dust 
ft CUT 


30 


20 


15 


10 


5 


Nevcnt 


1262 


1152 


1038 


849 


505 


P56 


12.7 


11.3 


10.1 


8.7 


6.7 


Acj) 


4.7 


4.3 


3.9 


3.2 


2.3 


D OUt 


10.0 


8.8 


7.8 


6.4 


4.7 


p t M>4.2 


2.0 


2.0 


2.0 


2.0 


2.0 




-0.0056 


-0.0074 


-0.0080 


-0.0055 


-0.0007 


(p t '-pr rt )/p t ' 


-0.0126 


-0.0135 


-0.0137 


-0.0120 


-0.0124 


{Pt 1 -Pt'MPt 1 


0.0054 


0.0042 


0.0039 


0.0049 


0.0098 


P t (0+77>5)/PJ 


-0.0006 


-0.0007 


-0.0001 


0.0022 


0.0083 


1 - cos{Acp) 


0.0060 


0.0050 


0.0040 


0.0027 


0.0014 


a{Db[ 7 ,J\) 


0.1207 


0.1012 


0.0897 


0.0743 


0.0620 


a{Db[^ 1 pari\) 


0.1442 


0.1212 


0.1083 


0.0937 


0.0806 


Entries 


7128 


6507 


5866 


4794 


2852 



Table 9: Selection 2. 90 < Pt 1 < 140 GeV/c {e^ et < 3%. 



r> ctust 
^tCUT 


30 


20 


[ 15 


10 


5 


Nevcnt 


797 


711 


632 


511 


288 


P56 


13.4 


11.6 


10.4 


8.9 


6.9 


Ac)) 


4.0 


3.4 


2.9 


2.4 


1.8 


l ) out 


10.9 


9.2 


8.0 


6.6 


4.8 


p t M>4.2 


1.9 


1.9 


1.9 


1.9 


1.9 


(PS-pr^/Pt 1 


-0.0100 


-0.0101 


-0.0092 


-0.0062 


-0.0018 


(P t J -Pt Part )/P t J 


-0.0160 


-0.0149 


-0.0137 


-0.0118 


-0.0105 


{Pt 1 -Pt^/Pt 1 


0.0045 


0.0036 


0.0034 


0.0047 


0.0077 


P t (0+??>5)/PJ 


0.0000 


0.0004 


0.0010 


0.0031 


0.0069 


1 - cos{Acp) 


0.0045 


0.0033 


0.0024 


0.0016 


0.0008 


a{Db[ 1 ,J\) 


0.0934 


0.0764 


0.0668 


0.0552 


0.0456 


a{Db[-f,part]) 


0.1145 


0.0956 


0.0872 


0.0763 


0.0624 


Entries 


17161 


15309 


13613 


11009 


6200 



Appendix 3 

pmin = 40 GeV/c 
P t lso1 < 2 GeV/c, < 5%, A0 = 17° (Selection 1) 



Table 1: Number of signal and background events (per Li n t = 300 pb 1 ). 



p cLust 
r t cut 

(GeV/c) 


P t %t (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


40000 


59000 


62000 


62000 


62000 


62000 


10 


50000 


96000 


112000 


115000 


115000 


115000 


15 


52000 


105000 


132000 


141000 


143000 


143000 


20 


53000 


107000 


139000 


153000 


158000 


159000 


30 


53000 


109000 


143000 


159000 


170000 


173000 



Table 2: S/B. 



p cLust 
"t cut 
(GeV/c) 


P tc Z (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


3.6± 0.3 


3.4± 0.2 


3.4± 0.2 


3.4± 0.2 


3.4± 0.2 


3.4± 0.2 


10 


3.3± 0.2 


2.9± 0.1 


2.8± 0.1 


2.8± 0.1 


2.8± 0.1 


2.8± 0.1 


15 


3.2± 0.2 


2.7± 0.1 


2.5± 0.1 


2.5± 0.1 


2.5± 0.1 


2.5± 0.1 


20 


3.1± 0.2 


0.0± 0.0 


2.4± 0.1 


2.4± 0.1 


2.4± 0.1 


2.4± 0.1 


30 


3.1± 0.2 


2.6± 0.1 


2.3± 0.1 


2.3± 0.1 


2.2± 0.1 


2.2± 0.1 


Table 3: (F), F = (Pj-P/ et )/PJ. 


p cLust 
r ' cut 

{GeV/c) 


P t ° c :t (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


0.012 


0.020 


0.023 


0.023 


0.023 


0.023 


10 


0.011 


0.024 


0.032 


0.034 


0.035 


0.035 


15 


0.011 


0.025 


0.035 


0.040 


0.042 


0.043 


20 


0.011 


0.025 


0.036 


0.041 


0.046 


0.046 


30 


0.011 


0.025 


0.035 


0.042 


0.047 


0.049 


Table 4: a(F), F = {P t ^-P t Jet ) / 'P{< . 


p cLust 
r t cut 

{GeV/c) 


P t ™t (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


0.053 


0.070 


0.074 


0.074 


0.075 


0.076 


10 


0.054 


0.080 


0.095 


0.099 


0.101 


0.102 


15 


0.055 


0.082 


0.104 


0.115 


0.121 


0.121 


20 


0.055 


0.083 


0.108 


0.123 


0.135 


0.137 


30 


0.055 


0.083 


0.109 


0.127 


0.150 


0.159 



pmin = 70 GeV/c 
P t lso1 < 2 GeV/c, < 5%, A0 = 17° (Selection 1) 



Tabic 5: Number of signal and background events (per Li„t = 300 pb 1 ). 



p clust 
r f cut 

{GeV/c) 


PtTt (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


2900 


4500 


4700 


4700 


4700 


4700 


10 


3600 


7100 


8500 


8900 


9000 


9000 


15 


3800 


7700 


10100 


11200 


11800 


11800 


20 


3800 


7900 


10600 


12300 


13600 


13700 


30 


3800 


8000 


10900 


12900 


15400 


16000 



Tabic 6: S/B. 



p dust 
r ( cut 

(GeV/c) 


P tc t (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


ll.lzb 1.1 


10.3± 0.8 


10.3± 0.8 


10. 2± 0.8 


10.1± 0.8 


10.0± 0.8 


10 


10.1± 0.9 


8.8± 0.5 


8.3± 0.4 


8.2± 0.4 


8.1± 0.4 


8.1± 0.4 


15 


9.8± 0.8 


8.2± 0.5 


7.4± 0.4 


7.1± 0.3 


6.8± 0.3 


6.8± 0.3 


20 


9.4± 0.8 


7.9± 0.4 


7.0± 0.3 


6.5± 0.3 


6.1± 0.2 


6.1± 0.2 


30 


9.3± 0.8 


7.6± 0.4 


6.6± 0.3 


6.0± 0.2 


5.4± 0.2 


5.2± 0.2 



Table 7: (F), F = (P^-P/^/P^. 



p clust 
"t cut 

(GeV/c) 


PtZi (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


0.008 


0.014 


0.016 


0.016 


0.016 


0.016 


10 


0.008 


0.015 


0.020 


0.023 


0.024 


0.024 


15 


0.008 


0.015 


0.022 


0.027 


0.031 


0.031 


20 


0.007 


0.014 


0.022 


0.028 


0.037 


0.039 


30 


0.007 


0.014 


0.022 


0.029 


0.043 


0.052 


Table 8: a(F), F = (Pj-P t Jet ) / 'P{< ■ 


p clust 
r t cut 

(GeV/c) 


P t %« (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


0.031 


0.042 


0.045 


0.046 


0.046 


0.046 


10 


0.032 


0.048 


0.058 


0.062 


0.064 


0.064 


15 


0.032 


0.049 


0.063 


0.072 


0.078 


0.078 


20 


0.032 


0.050 


0.065 


0.078 


0.089 


0.090 


30 


0.032 


0.050 


0.066 


0.080 


0.099 


0.102 



Pjnin = 1Q Q QeV / C 

P t iso1 < 2 GeV/c, < 5%, A0 = 17° (Selection 1) 



Tabic 9: Number of signal and background events (per L;„t = 300 pb 1 ). 



p cLust 
r t cut 
(GeV/c) 


P,SS (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


510 


820 


870 


870 


870 


870 


10 


630 


1270 


1560 


1630 


1650 


1650 


15 


650 


1380 


1830 


2050 


2150 


2150 


20 


660 


1410 


1930 


2260 


2520 


2560 


30 


670 


1430 


1970 


2370 


2870 


3060 


Table 10: S/P. 


p dust 
r f cut 

(GeV/c) 


P t ™' (GeF/c) 


5 


10 


15 


20 


30 


1000 


5 


29.5± 4.0 


26.5± 2.7 


25.3± 2.5 


24.9± 2.4 


24.9± 2.4 


24.5± 2.3 


10 


26.9± 3.1 


22.3± 1.7 


20.1± 1.3 


19. 0± 1.2 


18. 7± 1.2 


18. 7± 1.2 


15 


24.5± 2.7 


20.2± 1.4 


17.1± 1.0 


15. 7± 0.8 


14.9± 0.8 


14.9± 0.8 


20 


23.6± 2.5 


18.6± 1.2 


15. 5± 0.8 


13. 7± 0.7 


12. 3± 0.5 


12. 1± 0.5 


30 


23.1± 2.5 


18.3± 1.2 


14. 7± 0.8 


12. 7± 0.6 


10.8± 0.4 


10.0± 0.4 


Tabic 11: (F), F = (Pj-P t Jet )/P t ^. 


(GeV/c) 




5 


10 


15 


20 


30 


1000 


5 


0.007 


0.012 


0.013 


0.013 


0.013 


0.013 


10 


0.007 


0.012 


0.017 


0.018 


0.019 


0.019 


15 


0.007 


0.012 


0.017 


0.021 


0.024 


0.024 


20 


0.007 


0.012 


0.017 


0.022 


0.027 


0.029 


30 


0.007 


0.012 


0.017 


0.022 


0.030 


0.038 



Table 12: a(F), F = (P t ^-P t Jet ) / 'Pp . 



p dust 
r t cut 

{GeV/c) 


P t ™i (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


0.022 


0.031 


0.034 


0.034 


0.034 


0.034 


10 


0.023 


0.035 


0.042 


0.044 


0.045 


0.045 


15 


0.023 


0.035 


0.045 


0.052 


0.055 


0.055 


20 


0.023 


0.036 


0.046 


0.055 


0.061 


0.061 


30 


0.023 


0.036 


0.047 


0.057 


0.066 


0.067 



pmin = 4 Q Q eV / c 

P t lso1 < 2 GeV/c, < 5%, A0 = 17°, e jet < 3% (Selection 2) 



Tabic 13: Number of signal and background events (per L int = 300 pb 1 ). 



p dust 
r t cut 

{GeV/c) 


Ptlt (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


23000 


33000 


34000 


34000 


34000 


34000 


10 


27000 


47000 


53000 


54000 


54000 


54000 


15 


28000 


50000 


60000 


63000 


63000 


63000 


20 


28000 


51000 


62000 


66000 


68000 


68000 


30 


28000 


51000 


63000 


68000 


72000 


73000 



Tabic 14: S/B. 



p dust 
r t cut 

(GeV/c) 


P tc t (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


4.0± 0.4 


4.1± 0.4 


4.1± 0.4 


4.1± 0.4 


4.0± 0.4 


4.0± 0.4 


10 


3.9± 0.4 


3.7± 0.3 


3.6± 0.3 


3.6± 0.2 


3.6± 0.2 


3.6± 0.2 


15 


3.8± 0.4 


3.4± 0.2 


3.2± 0.2 


3.2± 0.2 


3.2± 0.2 


3.2± 0.2 


20 


3.8± 0.4 


0.0± 0.0 


3.2± 0.2 


3.1± 0.2 


3.0± 0.2 


3.1± 0.2 


30 


3.8± 0.4 


3.4± 0.2 


3.1± 0.2 


3.0± 0.2 


2.9± 0.2 


2.8± 0.2 


Tabic 15: (F), F= (P t x '-P t Jet ) / Pf 1 . 


p dust 
r < cut 

(GeV/c) 


P t ° c :t (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


0.006 


0.008 


0.008 


0.008 


0.008 


0.008 


10 


0.006 


0.005 


0.004 


0.003 


0.003 


0.003 


15 


0.006 


0.006 


0.003 


0.003 


0.001 


0.001 


20 


0.006 


0.006 


0.003 


0.000 


-0.002 


-0.003 


30 


0.006 


0.006 


0.003 


0.001 


-0.005 


-0.006 


Tabic 16: a(F), F = (P^-P t Jet )/P t \ 


p dust 
r < cut 
{GeV/c) 


P t %* (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


0.050 


0.066 


0.069 


0.069 


0.069 


0.069 


10 


0.052 


0.074 


0.086 


0.089 


0.090 


0.090 


15 


0.051 


0.075 


0.095 


0.102 


0.107 


0.107 


20 


0.052 


0.075 


0.097 


0.109 


0.120 


0.123 


30 


0.052 


0.075 


0.098 


0.113 


0.136 


0.147 



pmin = 70 GeV/c 
P t lso1 < 2 GeV/c, < 5%, A0 = 17°, e jet < 3% (Selection 2) 



Tabic 17: Number of signal and background events (per Li n t = 300 pb 1 ). 



p cLust 
r f cut 
(GeV/c) 


P t °£ (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


2300 


3400 


3600 


3600 


3600 


3600 


10 


2800 


5000 


5800 


6000 


6000 


6000 


15 


2900 


5300 


6700 


7200 


7400 


7400 


20 


2900 


5400 


6900 


7700 


8200 


8300 


30 


2900 


5500 


7000 


8000 


9000 


9200 



Table 18: S/B. 



p clust 
r ( cut 
(GeV/c) 


P t <Z (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


11.8± 1.3 


11.6± 1.1 


11.4± 1.0 


11. 3± 1.0 


11. 3± 1.0 


ll.lzfc 1.0 


10 


ll.Ozfc 1.1 


10.1± 0.8 


9.5± 0.6 


9.3± 0.6 


9.3± 0.6 


9.3± 0.6 


15 


10.9± 1.1 


9.6± 0.7 


9.0± 0.6 


8.6± 0.5 


8.4± 0.5 


8.4± 0.5 


20 


10.6± 1.1 


9.3± 0.6 


8.5± 0.5 


8.1± 0.5 


7.5± 0.4 


7.4± 0.4 


30 


10.5± 1.0 


9.1± 0.6 


8.2± 0.5 


7.6± 0.4 


6.8± 0.3 


6.6± 0.3 



Table 19: (F), F = {P{ l -P t Jet )/P{i . 



p clust 
"t cut 

(GeV/c) 


P t %i (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


0.006 


0.009 


0.009 


0.009 


0.009 


0.009 


10 


0.006 


0.007 


0.008 


0.009 


0.009 


0.009 


15 


0.005 


0.007 


0.009 


0.009 


0.009 


0.009 


20 


0.005 


0.006 


0.008 


0.009 


0.010 


0.010 


30 


0.005 


0.006 


0.008 


0.009 


0.010 


0.010 


Table 20: a(F), F = (P t ^-P t Jet ) / ' P t \ 


p clust 
r t cut 

(GeV/c) 


P t ™t (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


0.031 


0.040 


0.043 


0.043 


0.043 


0.043 


10 


0.031 


0.046 


0.054 


0.057 


0.058 


0.058 


15 


0.031 


0.047 


0.059 


0.066 


0.069 


0.069 


20 


0.031 


0.047 


0.060 


0.071 


0.078 


0.078 


30 


0.032 


0.047 


0.061 


0.073 


0.086 


0.088 



pmin = 10Q Q eV / c 

P t lso1 < 2 GeV/c, < 5%, A0 = 17°, e jet < 3% (Selection 2) 



Table 21: Number of signal and background events (per Lj n ( = 300 pb 1 ). 



p cLust 
(GeV/c) 


P t ™J (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


460 


720 


760 


760 


760 


760 


10 


560 


1060 


1250 


1300 


1300 


1300 


15 


580 


1130 


1440 


1570 


1620 


1620 


20 


580 


1150 


1490 


1700 


1830 


1840 


30 


580 


1160 


1520 


1750 


2020 


2090 



Table 22: S/B. 



p clust 
r f cut 
(GeV/c) 


P t 2t (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


31.9± 4.7 


27.6± 3.0 


26.6± 2.8 


26. 0± 2.7 


26.0± 2.7 


25. 5± 2.7 


10 


31. 1± 4.1 


23.9± 2.2 


22.4± 1.7 


21. 7± 1.6 


21. 5± 1.6 


21.4± 1.6 


15 


29.5± 3.7 


22.6± 1.8 


19.1± 1.3 


18. 2± 1.2 


17.5± 1.1 


17.5± 1.1 


20 


28.7± 3.6 


21. 5± 1.5 


17.5± 1.1 


15. 9± 0.9 


14.7± 0.8 


14.6± 0.8 


30 


28.1± 3.5 


20.9± 1.6 


16.6± 1.0 


14.7± 0.8 


12. 9± 0.6 


12.3± 0.6 



Table 23: (F), F = {P{<-P t Jet )/P{i '. 



(GeV/c) 




5 


10 


15 


20 


30 


1000 


5 


0.006 


0.009 


0.010 


0.010 


0.010 


0.010 


10 


0.006 


0.009 


0.010 


0.011 


0.011 


0.011 


15 


0.006 


0.008 


0.010 


0.011 


0.011 


0.011 


20 


0.006 


0.008 


0.010 


0.010 


0.011 


0.011 


30 


0.006 


0.008 


0.009 


0.010 


0.011 


0.011 


Table 24: a(F), F = (P^-Pt Jet )jP{ 1 . 


tj clust 
"t cut 
{GeV/c) 


P t ° c :t (GeV/c) 


5 


10 


15 


20 


30 


1000 


5 


0.022 


0.030 


0.032 


0.032 


0.033 


0.033 


10 


0.023 


0.033 


0.040 


0.042 


0.042 


0.042 


15 


0.023 


0.034 


0.043 


0.048 


0.050 


0.050 


20 


0.023 


0.035 


0.044 


0.051 


0.056 


0.056 


30 


0.023 


0.035 


0.044 


0.053 


0.060 


0.060 



